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ABSTRACT 

This  thesis  presents  some  general  equations  describing  the 
particle  density  distributions  vithin  a  plasma.  Solutions  to 
these  equations  are  discussed  for  various  geometries,  vith  particular 
emphasis  given  to  the  effects  of  diffusion  losses  at  probes  of  dimen- 
sions comparable  to  or  greater  than  the  mean  free  path  and  to  magnetic 
field  end  effects.  For  many  of  the  geometries  discussed,  including 
one  probe  geometry,  the  effects  of  a  superimposed  magnetic  field  is 
discussed. 


INTRODUCTION 

One  of  the  complex  problems  in  the  field  of  plasma  physics  is 
that  of  predicting  the  effects  of  surfaces  vith  vhich  the  plasma  is  in 
contact  upon  the  particle  density  distribution  functions  and  related 
phenomena.  No  plasma  is  infinite,  and  consequently  a  plasma  must  be 
confined  by  vails,  vhich  may  be  either  conductors  or  insulators. 
Further,  in  order  to  make  measurements  vithin  a  plasma,  probes  may 
be  inserted.  Vhat  is  their  effect  on  the  particle  density  distri- 
bution and  related  phenomena  such  as  potential  distribution  and 
energy  dissipation?  that  are  the  particle  currents  to  various  surfaces? 
The  latter  are  particularly  important  because  they  determine  the 
heat  liberated  at  the  surface,  vhich  in  turn  is  an  energy  loss  vhich 
may  cause  very  high  temperatures,  increased  local  ionization,  and 
melting  of  materials  in  extreme  cases.  This  thesis  discusses  these 
effects  vhen  the  surfaces  are  insulators. 

In  this  thesis  some  general  equations  describing  a  plasma  are 
presented  vith  their  corresponding  restrictions.  Solutions  to 
these  equations  are  discussed  for  various  geometries,  vith  particular 
emphasis  given  to  the  effects  of  probes  and  to  magnetic  field  end 
effects  on  the  steady-state  density  distributions  and  on  the  energy 
losses  at  the  surfaces  involved.  For  many  of  the  geometries  dis- 
cussed, including  one  probe  geometry,  the  effects  of  a  superimposed 
magnetic  field  is  studied. 


II 

GENERAL  RELATIONS 
In  this  section,  general  equations  describing  a  plasma  are 
developed  and  discussed. 
A.  Conservation  Equations 

The  Boltzmann  Transport  Equation 


is  used  as  the  starting  point  for  the  development  of  the  macroscopic 
relations  governing  the  plasma.  By  ascribing  appropriate  values  to  X 
and  >'  [I]   [2],  we  obtain  the  continuity  equation; 

f*  4-V-frw^  =  riC^-Pj  (2-2) 

(in  which  £[•  is  the  mean  ionizing  frequency  and  ;^  the  mean  recombination 
frequency)  ,  and  the  equation  for  conservation  or  momentum: 

|K0  t^f  -«(J  +  ^v^  =  -m^  (2-3) 

in  which  the  right  hand  term  is  really  a  definition  of  the  mean 
collision  frequency,   Vc   ,   the  pressure    *p     is  a  symmetric  tensor 

pj.Tnm^fy  (2-3a) 


and   in  which 


and 


a=4rE  C2-3b) 


m 


tJw-  -^S.  (2-3c) 

E,  B,  q,  m,   and  vfj  pre,   respectively,   the  electric  field   intensity, 
the  magnetic  field  intensity,  particle  charge,  particle  mass,   and 
drift  velocity. 


The  continuity  equation  (2-2)  can  be  further  simplified.     We 
shall  consider  a  tvo  component  system  of  electrons  and  singly-charged 
ions   in  which  volume  ionization  by  electrons  is  essential.     Recombina- 
tion takes  place  predominately  at  the  vails,  with  volume  recombination 
being  negligible,     VJe  have,  then, 

vhere  9  is  the  net   ionizing  frequency  and  the  particle  current  density 
is 

— »■ 

Pt  '-  Ykfa.  (2-La) 

In  the  above  equations,  and  in  the  remainder  of  this  thesis,  the 
subscript  "+"  refers  to  singly-charged  positive  ions,  and  the  sub- 
script "-"  refers  to  electrons.  Under  those  conditions  for  which 
the  rates  of  decay  of  ions  and  electrons  are  essentially  equal,  so 
that 

we  obtain  for  the  continuity  equation 

V-fl*  V^  =P77.  -fg  (2-6) 

In  the  momentum  equation  (2-3),  the  pressure  term  is  generally 
a  difficult  one  to  evaluate  because  "p  is  not  known,  and  the 
dependence  on  its  divergence  makes  good  approximations  difficult. 
However,  for  the  simple  case  in  which  the  distribution  of  random 
velocities  is  essentially  isotropic,  so  that 

V^/'^O  L/j,  (2-7a) 

and 

V^1  *  Vj"   »  V37,  (2-7b) 
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the  pressure  tensor  reduces  to 

p*  -3-  -n*T  =   -jri  (2.8) 

where  B  and  yU  are  the  diffusion  and  nobility  coefficients  respectively, 
The  last  term  is  derived  from  the  preceeding  ones  by  the  use  of 
Einstein's  relations 

This  generalization  (equation  (2-8))  is  valid  only  for  a  plasma  which 
is  characterized  by  a  small  drift  velocity, 

|Va\«J^.  (2-10) 

The  last  term  in  equation  (2-3)  is  frequently  written 

.y,ujA   =  -  ^5.  (2-11) 

which  defines  the  mobility,  la   .  In  the  case  of  ion-neutral  particle 
interactions,  ^  is  a  constant  because  the  forces  are  inverse  fifth 
power  polarization  forces.  In  all  other  cases  we  must  introduce  an 
appropriate  average  value,  V   ,  which  we  w£ 11  consider  independent 

of  v.. 

d 

One  important  consequence  of  the  above  assumptions  is  that  the 
ion  and  electron  mobilities  and  diffusion  coefficients  are  scalar 
constants. 

Substitution  of  equations  (2-8)  and  (2-11)  into  the  momentum 
equation  (2-3),  and  rearranging  terms,  yields,  for  the  positive 
ions  of  mass  Ms 

and  for  the  electrons  (mass  m); 


For  stationary  (steady-state)  situations  in  the  plasma,  ve  find 
that  these  equations  reduce  to 


B.     Current  Density,  Field  and  Potential  Equations 
From  Poisson's  equation  we  have 

and  from  Maxwell's  equations  for  the  stationary  case  2 

V*E=0;  (2.l6) 

These  three  equations,  with  the  continuity  equation  (2-6)  and  the 
simplified  momentum  conservation  equation  (2-lIt)  express  the 
governing  conditions  for  a  plasma  vhich  complies  with  the  assumptions 
ennunerated  to  date.  These  will  be  solved  for  current  equations  and 
for  general  field  and  potential  relations. 

From  equation  (2=6)  it  is  readily  deduced  that 

Fl   .P.  -rfVxP)  (2=18) 

— ^       — * 

in  vhich  (V  x  F  )  is  a  vector  vhich  we  shall  now  determine 
Let  us  define  three  orthogonal  directions  in  the  plasma,  one  parallel 
to  the  magnetic  field,  designated  by  the  subscript  "  ||  ",  one 
perpendicular  to  the  field  and  to  any  surfaces  parallel  to  the  field, 
"  J^",  and  the  third  transverse  to  the  field  and  to  the  perpendicular 


direction  "T" .  In  cylindrical  coordinates,  with  a  field  parallel  to 
the  z-axis  and  cylindrical  boundaries  parallel  to  the  field  these 
would  be  the  z-j  r-5  and  if?-    directions  respectively.  The  ion  and 
electron  particle  currents  at  the  sheath  edge  are  equal  at  all  points 
for  a  steady-state  plasma  because  they  must  arrive  in  equal  numbers 
to  recombine.  For  a  geometry  in  which  the  insulator  surfaces  coicide 
with  coordinate  surfaces  (the  only  geometries  for  vhich  we  will  be 
able  to  obtain  exact  solutions  for  the  particle  density  distributions), 

rv  -  r\  =  vL  (2-1*0 

and 

r^-r.,,.^  (2-1*0 

at  the  sheath  edges.  We  now  impose  the  condition  that 

VTvu   *^r\-^£T=o„  (2-20) 

This  requirement  is  met  either  by  a  plasma  vhich  has  a  rotational 
axis  of  symmetry  parallel  to  the  magnetic  field  or  one  in  which  the 
field  lines  coincide  with  one  set  of  orthogonal  coordinate  lines  and 
in  vhich  the  resulting  distribution  is  a  function  of  only  one  coordinate. 
The  latter  condition  would  be  met  in  the  vicinity  of  a  plane  surface 
parallel  to  the  magnetic  field  lines  in  an  "infinite"  plasma.  Here, 
for  a  y-z  plane  and  a  field  parallel  to  the  z-axis,  the  perpendicular 
direction  would  be  the  x-direction  and  the  transverse  direction  the 
y~direction.  If  these  symmetry  conditions  are  imposed,  it  is  readily 
deduced  that  equations  (2-19)  apply  throughout  the  plasma  if  the 
walls  are  insulators.  If,  however,  the  end  surfaces  are  conductors 
(i.e.,  electrodes),  equation  (2-19b)  does  not  apply.  The  parallel 
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current  is  then  a  conduction  current,  not  a  diffusion  current,  and  is 
not  treated  by  the  developments  to  follow.  This  current  does  not 
affect  the  perpendicular  currents,  however,  and  such  a  plasma  will 
be  treated  as  an  infinite  plasma  in  the  longitudinal  direction,  with 

V»n  -  O.  (2-21) 

From  equation  ( 2— lit)  and  the  symmetry  conditions  (2-20),  it  is 
immediately  seen   that  the  transverse  component  equation  is: 

rVy-  =  +/u+pj3.  (2-22) 

Fropi  this  equation,  and  equations  (2-19),  we  can  see  that     (v*t"  / 
of  equation  (2-18)  has  only  a  transverse  component: 

(VxF)   =fVxF)T=    -(,,#)  Cx  (2-23) 

The  parallel  and  perpendicular  component  equations  derived  from 
equation  (2-lh)  are 

Pi,,  --  */>!«»*£„ -Dt^n*,  (J-2« 

and 

V+i   +/JiStT&  =  tytrttEj.   -Dt^n*.  (2-25) 

A  simple  solution  is  possible  if  the  assumption  of  proportionality 
holds: 

%*%-=^  (M6> 

By  eliminating  Em  from  equation  (2-2h),  we  obtain 

vhere  Dm.  is  a  quasi -ambipolar  diffusion  coefficient  for  diffusion 
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parallel  to  a  magnetic  field  and  Ds  is  the  ambipolar  diffusion 
coefficient  which  is  commonly  developed  for  a  plasma  vhen  B  =  0  by- 
infer  ing  congruence, 

from  the  current  density  divergence  equality  equation  (2-6)  and  from 
the  proportionality  condition  (2-26)  [2~]  [3D  •  Ds  *s  related  to 
Schottky»s  ambipolar  diffusion  coefficient, 


by  the  equations 


U^  -    U<a.   

fax*  +  M'*'i 


(2-30) 


-DcO^^).  (2.30a) 

In  equation  (2-30a)  vre  have  employed  the  following  relations  for 
charge  density  and  conductivity 

p±bl+-7l')e  (2-3  ia) 

and 

CT=  fair +/**+&.  (2°31b) 

By  substituting  equation  (2-22)  into  equation  (2-25),  we  obtain 

fVj.  +  /utB^r*i.  -  */****  E^-D^ti*.  (2-32) 

By  the  use  of  equation  (2-19a),  E,  may  be  eliminated  from  the 
electron  equation,  to  yield 


By  simplifying  and  solving  for  \A    ,  we  obtain 


"  D5\7±77- 


n  "    r  ~"^'  -,.  fl "  -^V^y)-,     (2-3I4) 


/  + 


,  4,  ^^--/^) 


/^-  Bz 


wnere  Dm,  is  a  quasi -ambipolar  diffusion  coefficient  for  diffusion 
perpendicular  to  a  magnetic  field.  Thus  diffusion  perpendicular 
to  the  magnetic  field  is  reduced  by  an  inverse  dependence  on  the 
square  of  the  magnetic  field  intensity. 

If  there  is  no  superimposed  electric  field,  so  that  E  -   Es, 
vhere  Es  is  the  space-charge  field,  ESj_  may  be  obtained  by  substituting 
equation  ( 2-3^)  into  (2-32); 

Es,  =  r~ ^  (2-35) 

J-         /A.  >1 

Similarly, 

F<=  -  D^rD"  %iH  -   ^s-D-  %n  (2-36) 

The  transverse  spac^  charge  field  is  zero  as  a  consequence  of  our 
assumptions  that  there  is  no  transverse  density  gradient,. 
For  a  quasi -neutral  plasma,  in  which 

it  follows  from  equation  (2-30)  that  Ds  «  Da,  and  that 

In  the  absence  of  a  magnetic  field,  equation  (2-35)  when  simplified 
and  added  to  equation  (2-36),  yields  the  single  equation  for  the 
space-charge  field 

f    =  Ps-D-  %_  D^-q    V„ 

Ls    /A-        0       iv»>1-  m*+m-      ^  '       (2-39) 


From  the  field  equations,  we  obtain  the  following  equation  for 
the  space  charge  potential  for  a  plasma  with  only  a  perpendicular 
density  gradients 


V,  = 


D-~(l+>u.8 


(2-10) 


In  the  absence  of  a  magnetic  field,  equation  (2-39)  yields 


y.  ^.a#-\ ->.   V  **^(£)  (2-u: 

The  zero  for  potential  in  both  of  these  equations  is  taken  as  the 
"center"  of  the  plasmas  that  point  (or  region)  where  the  density 
gradient  is  zero. 


C.  The  Diffusion  Tensor 

The  current  density  equations  (2-22),  (2-27)  and  (2»3ii)  can  be 
combined  into  one  equation 

[\    =  -D^±Vn_;  (2-Ji2) 

which,   in  rectangular  coordinates,    is 


D 


mx      Dm,T     O  \/  %">- 


--    D 


rv,aT 


t> 


-r/1, 


O     O   Dm, 
and  in  cylindrical  coordinates  is 


(2-!;3) 


<2-!ii) 


The  components  of  the  diffusion  tensor  are 


; 


EWT  -  -/U-Bl^'.  (2-l£d) 

We  shall  now  examine  each  of  these  components.  Equation 
(2»hf>a)  reflects  the  fact  that  the  magnetic  field  does  not  affect 
diffusion  in  the  direction  parallel  to  itself.  Further,  we  have 
shown  that 


"^MM  ~  ^5  f>->o     >  ^  •  (2~l;6a) 

ct  >u  *  ^  - 

A  His  [2]  has  also  shown  that 

Qs >D-,  (2-Wb) 

where  ^  is  the  conductivity  at  the  center  of  the  plasma.  The  quasi*- 
neutral  limit  (3*-ii6a)  is  applicable  to  a  plasma  in  which  the 
mobilities  (and  consequently  the  conductivity)  are  so  great  that 
negligible  space  charge  can  develop.  The  above  limits  are  the  low 
and  high  limits  for  the  ambipolar  diffusion  coefficient. 

The  perpendicular  component  is  more  complex.  First,  as  the 
magnetic  field  approaches  zero, 

H    ^  n  ( 2-Ii7a) 
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As  B  becomes  very  large,  Dm,  approaches  zero  as  "g* 

TN    ^  _E^ >£>  (  2-U7b) 

^•\_  g->«^  //^-  6*-   B-*«»    - 
From  equation  (2-4i6a)  it  is  seen  that 

and  from  (2-60a)  that 

D^^r  ^.^v^-s-  ,  (2-h7d) 

which  can  be  non-zero  only  if    j00-  0.     As  shown  in  equation  (2-5>!?)> 
for  a  long,  narrow  plasma  this  can  be  true  only  if 

Dr»,7 >  ,    ^~.a>  (^Te) 

From  these  equations  it  is  apparent  that  for  the  parallel  and 

perpendicular  diffusion  tensor  components  in  a  magnetic  field  the  high 

and  low  limits  are  the  corresponding  electron  diffusion  values  or  the 

quasi-neutral  ambipolar  values  respectively. 

For  the  transverse  components,  we  obtain 


+■  jut.  D*B 


& 

— *  Cj 

XV- 

B 

—7  CO 

->- 

(2-ii8a) 


(2-U8b) 


^O  (2-lrfc) 


B 


^/U±T^B    ^.^  »  O,  (2-li8d) 
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The  maximum  transverse  internal  current  will  flow  when  Dnu-p  are  a 
maximum.  For  the  quasi-neutral  plasma,  equation  (2-l|8b),  this 
occurs  ■when  the  magnetic  field  has  a  magnitude  given  by 

For  helium,  using  data  from  reference  \Jk\   ,  this  is  approximately  0.12 
weber/m  .  The  magnitude  of  the  current  density  resulting  from  this 
field  can  be  determined  from  equation  (2-22) 

vhere  £  is  a  unit  vector  in  the  transverse  direction.  The  magnetic 
field  generated  by  this  current  vi 11  oppose  the  applied  field,  and 
accounts  for  the  observed  diamagnetism  of  the  plasma.  This  is  the 
term  which  satisfies  equation  (2-17),  and  it  must  be  small  compared  to 
B  to  justify  neglecting  its  effects  on  the  other  plasma  equations. 

In  the  quasi-neutral  plasma,  with  ?74  ^>f__  ,  the  diffusion  tensor 
reduces  to 


Itf4t/U-B 


1  .*- 

1+fkp-B* 

O 

D.  The  Continuity  Equation 

By  substituting  equation  (2-l|2)  into  the  continuity  equation  (2-6), 
we  obtain 

§£  - V-(DmtVnJ)  =  \)r\.  (2-52) 
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a  second  order  partial  differential  equation.  For  the  steady  state, 
this  reduces  to 

whereas  for  a  decaying  plasma  (for  which  \i  ■  0),  it  becomes  the 
diffusion  equation 

||  _  f^& .)  =  o  (2.51;) 

In  order  to  simplify  equation  ( 2-510  >  we  must  examine  its  com- 
ponents. The  components  of  the  diffusion  tensor,  equations  (2-li5), 
are  constants  because  of  our  assumptions  and  consequently  the 
divergence  of  the  diffusion  tensor  is  zero.  In  the  absence  of  a 
magnetic  field,  D^+is  the  scalar  D  ,  and  equation  (2-52)  becomes  the 
familiar  continuity  equations 

For  the  steady  state,  we  can  further  simplify  this  by  employing  the 
assumption  of  proportionality,  equation  (2-26),  to  obtain 

n  v$ 

In  the  presence  of  a  magnetic  field,  we  shall  limit  ourselves  to 
orthogonal  geometries  in  vhich  equation  (2-53)  can  be  separated  into 


.•^viV).  4-  D*hVnn~  +  »n„*0  (2-5?) 

In  equation  (2-57)  we  have  made  explicit  use  of  the  requirement  that 
the  transverse  density  gradient  be  zero.  By  employing  the  assump- 
tion of  proportionality  (equation  (2-26))  and  a  theorem  from  vector 

analysis,  we  see  that 

lh 


: 


vith  a  similar  equation  for  the  perpendicular  components.  Thus  ve 
can  divide  equation  (2-57)  by  n_  and  expand  it  to  obtain 

If  ue  then  define  a  vector  u"  such  that 


— »  ■  ■  - 

ve  obtain 

vhich  is  a  first  order  partial  differential  equation  in  u.  In  the 
absence  of  a  magnetic  field,  D  is  a  scalar  (D  ),  and  equation  (2-56) 
can  be  similarly  expressed  (in  any  coordinate  system)  as 


V<U     +  U*  +£=0.  (2-62) 


The  significance  of  the  transformation  of  this  equation  from  a 
second  order  differential  equation  in  particle  density,  n,,  to  a 
first  order  differential  equation  in  u  will  become  apparent  in  the 
discussion  of  boundary  conditions  in  a  later  section  • 
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Eo  Boundary  Conditions 
1.  Derivation 

We  shall  now  develop  an  equation  for  the  boundary  condition 
and  evaluate  it  for  a  few  simplified  cases.  We  consider  a  quasi- 
neutral  stationary  plasma  in  contact  with  a  plane  wall.  The  mag- 
netic field,  if  any,  is  parallel  to  the  wall. 

There  are  two  conditions  which  we  must  fulfill.  The  first 
is  that  the  electron  and  ion  currents  to  the  insulator  surface 
must  be  equalj  the  other  is  that  the  current  must  be  continuous  at 
the  limit  of  the  diffusion  region,  which  can  be  at  the  edge  of  the 
sheath  or  at  the  edge  of  the  free-fall  zone. 

Ye   have  to  distinguish  different  sheath  theories  accord- 
ing to  whether  D-,  >  D+J_  ,  D-±  <  D+J_  or  B-±  ^  D^^  ,  where 


Du   -    T^  (2"66) 


a.     The  case  D-j_>>    D^j_       ,    ion  current  inertia-limited. 

Let  us  first  examine  the  problem  in  which  D_j_  >>      D+J_ 
In  this  case  a  negative  wall  charge  builds  up  until  the  electron 
current  is  reduced  to  equality  with  the  positive   ion  current. 
The  electrons  will  be  treated  as  having  a  Maxwell ian  distribution 

r\.  «r?ee  kT"  (2„67) 
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where  the  potential  V  at  the  sheath  edge   is  taken  as  zero,   -e   is   the 
charge  of  the  electron,,  and  the  subscript  "o"  means  evaluation  at  the 
sheath  edge.     This  notation  vill   be  used  throughout   this  section. 

Ye  find  that   two  possibilities   exist  for  the   ion  motion  through 
the  sheathi    (1)    their  motion  may  be   inertia-limited,   or  (2)    their 
motion  may  be  mobility-limited.      In  the   inertia-limited  case,    the 
sheath   is  a  free-fall  region  and  the  current   in  the  sheath   is   that 
given  by  kinetic  theory?  * 

r   =  ££s  (2.68) 

We  must  next  evaluate  the  diffusion  current  to  the  sheath 
edge,  which  is 

P   * ^ ^Vn  (2-69) 

If  the  electrons  and  ions  are  in  equilibrium  within  the  body  of  the 
plasma 

Da.  *  2 IV     (XW),  (2-69a) 

In  order  that  the  number  of  particles  varried  to  the  sheath  by 
the  above  diffusion  current  equal  the  number  carried  to  the  vail  as 

expressed  by  equation  (2-68),  the  normal  component  of  the  diffusion 
current  must  equal  the  sheath  current.  Therefore,  we  have  the 
equation 

noK   =  -tO^Vnn.  (2-70) 

From   this  equation  we  see  that 


a  result  obtained  by  Ecker  [3~],  where,   for   this  case 
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oL  =   jfO +/*+/*-&*).  (2-70a) 

It  will  be  seen  that  the  boundary  condition  can  alvays  be  ex- 
pressed in  the  form  (2-71),  with  the  value  of  oC  varying  with  the 
sheath  characteristics.  Let  us  investigate  the  potential  function  V 
for  this  condition.   If  there  is  no  magnetic  field,  ve  see  from  the 
conservation  of  energy  that 

M  V  .  Mv*   eV  (2^72) 

where  M  is  the  ionic  mass.  From  this  ve   can  derive  the  following 
equation  for  the  ion  current  density; 


Let  us  nov  employ  a  coordinate  system  in  which  the  x-direction 
is  normal  to  the  sheath  and  positive  towards  the  wall,  and  whose 
origin  is  at  the  sheath  edge.  Solving  equations  (2-67)  and  (2-73) 
for  n_  and  n+  respectively,  and  substituting  in  Poisson's  equation 
(2-15),  ve  find 


(2-73) 


Multiplying  by  3~  and  integrating  yields:  t 

where  we  have  set   -r-J   ■  -E^ and  Vj  =  0.  This  result  has  been 
reported  by  Ecker[5J«  For  the  sheath  to  be  stable,  -<---£  must 
increase  with  x.  Therefore  the  quantity  in  the  brackets  of  equation 
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(2-75)  nust  be  positive.  Expanding  this  quantity  in  series,  we   find 


eV      ier    .   *tw  .Ley       iii\  >n  7AN 


feT- 

from  vhich  ve   see  that 


-r     >  -i  T 
/+     -  a  i-  (2-77) 

is  the  condition  for  the  sheath  to  be  stable.  Bohm  (jQ  has  calculated 

the  function  V  for  two  values  of  T+  by  numerical  integration. 

If  we  neglect  the  electron  density  in  Poisson's  equation,  and 

integrate,  we  obtain  the  Langmuir  relation  C83s 

V=-(^(ft*f3  (2-78, 

b.  The  case  D-j_.»  D-t-j,  ,  ion  current  mobility- United 

We  shall  next  examine  the  case  in  which  the  ion  current 

through  the  sheath  is  mobility  limited;  i.e.,  the  ions  move  in  the 

field  of  the  negative  wall  charge  according  to  their  mobility.  Nov 

the  electron  current  density  striking  the  wall  is  obtained  from  kinetic 

theory  and  the  electron  density  relation  (2-67) s 

.  eV(d.) 
-\  -7  k-T- 


*  /cL  '  4-  >  (2-79) 

where  V(d_)  is  evaluated  at  the  edge  of  the  electron  free-flight  zone 
(suitably  corrected  to  account  for  curvature  of  the  electron  paths  in 
the  magnetic  field).  Again  equating  this  current  with  the  diffusion 
current  brought  to  the  sheath,  we  obtain? 
eVCi.) 
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In  order  to  further  simplify  this  relation  for  the  case  in  which  the 
electrons  and  ions  are  in  equilibrium,  ve  use  the  relation 

mvj1  =  M7+*  (2-81) 

from  which  ve  see   that 

^  -  ~  n,  M  e   Kr-  (i  ^/,-S1) .  (2-82) 

Thus  we  obtain  for   the  parameter    ^      : 

In  order  to  evaluate   this  expression,    it  is  necessary  to   knov  the 
potential  function  V(x).     The   ion  current  density   in  the  sheath 
(neglecting  the  magnetic  field)    iss 

y  -  /V'e E°  =  M+i/i+e Z*  -  -  /**.>!* .e  ^  ?  ( 2-81) 

where 

,,  __Z^—  (2-81*) 

Again  employing  Poisson's  equation,  and  solving  equations  (2-8ii)   and 
(2-67)   for  n+  and  n-  respectively,   we  obtain: 

At  the  sheath  edge  (where  x  =  0,  V  =  0,   and     £p    =  +E  ),    ^p.     -  0, 
a  result  consistent  with  the  assumption  of  quasi-neutrality  at  the 
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d2V 

e  f     sY 

+ 

4- 

sheath  boundary.  Equation  (2-85b)  can  be  integrated  once,  and  yields 
a  first  order  differential  equation  for  V; 

Further  integration,  to  obtain  V(x),  must  be  performed  numerically. 
However,  if  we  again  neglect  the  electron  density,  equation  (2-86) 
can  be  integrated  to  yield? 


J  e0ju+  (2-87) 

c.  The  case  D^i^  D+l 

An  approximation  for  the  case  D-x  »  d+j_   can  be  easily 
obtained  for  inertia-limited  ion  current.  By  neglecting   ^Mf2(52 
as  compared  to  1  in  the  equality 


we  obtain 


(2-88) 


!2~    *    I  +  {A*t ",  (2-89) 


from  which  it  is  seen  that 


f-(vt-*}     -'^-B%  (2.90) 

Employing  this  relation  in  equation  (2-?0a)  yields 


If  we  neglect  the  last  term  in  the  parenthesis  as  compared  with  one, 
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and  employ  Einstein's  relation 
vz   obtain 


(2-92) 


^  -  g(l  +  tX  (2-93) 

do  The  case  D+j^D..^ 

For  the  plasma  in  which  D+j>>  D«^  t   a  positive  wall  charge 
must  build  up  until  the  ion  and  electron  currents  are  equal.  The  ions, 
assuming  Maxvellian  distribution,  vill  have  a  density  within  the  sheath 
given  by 

m^-yioe  *T%  (2-910 

where  the  ionic  charge  is  +e„ 

The  ion  current  to  the  vail  is  obtained  from  kinetic  theory  as 

where  the  evaluation  of  V(d+)  is  at  an  ionic  mean-free  path,  d+ , 


which  has  been  adjusted  to  include  the  effects  of  the  magnetic  field. 
Equating  this  to  the  diffusion  current  to  the  sheath  yields  the 
boundary  condition  (2-71)  in  vhich 

In  order  to  evaluate  the  function  V,  it  is  necessary  to  employ 
the  equation  for  the  electron  current  (neglecting  the  effects  of 
the  magnetic  field)? 
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By  the  sane  procedure  used  for  the  mobility- limited  ions  vhen 
D->>D+J_,  we  obtain 

(Iff '%[-&<■  W*-k*&        («e, 

2 .  Evaluation  of  ot 

The  over-all  problem  of  evaluating  ot  of  tbe  boundary  condi- 
tion equation  (2-71)  is  a  much  more  complicated  problem  than  indica- 
ted by  the  derivations  of  tbe  previous  section.  The  influence  of  the 
transition  region,  the  effect  of  non-planar  geometries,  deviation 
from  thermal  equilibrium  within  the  plasma  proper  (and  in  the  transi- 
tion region),  and  secondary  emission  are  just  a  few  of  the  factors 
which  can  affect  this  quantity.  The  transition  region,  which  joins 
in  a  continuous  fashion  the  sheath  and  the  main  body  of  the  plasma, 
presents  particularly  difficult  problems  which  are  beyond  the  scope 
of  this  thesis.  In  the  computation  of  qC   for  a  surface  whose  dimen- 
sions are  comparable  to  or  less  than  the  mean  free  path,  it  is 
especially  important  to  compute  the  fraction  of  particles  brought 
to  the  sheath  edge  which  would  miss  the  surface  in  their  flight 
through  the  sheath  and  would  pass  back  into  the  plasma.  Thus,  for 
the  limit  where  the  surface  is  very  much  smaller  than  tbe  mean  free 
path,  as  for  a  Langmuir  probe,  ot  can  become  vanisbingly  small. 

3.  Discussion 

An  understanding  of  the  processes  occuring  within  the  sheath 
has  been  shown  to  be  necessary  for  the  evaluation  of  the  boundary 
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condition  equation  (2-71).     We  have  seen  that  the  requirement  that 
the  current  be  continuous  results   in  this  condition  t-hich  must  be 

met  at   the  surface  of  the  ion  free-fall  zone  (the  region  within  one 
ionic  mean  free  path,     X+  ,  of   the  vail)  or  at   the  surface  of  the 

sheath,   whichever   is  farther  from   the  vail.     As  the  sheath  thickness 
and      Xir    become  very  much  smaller    than  other  dimensions,   equation 
(2-71)  may  be  replaced  by   its   limit  for  vanishing       )y+  ; 

*)^  =°-  (2-99) 

This  is   the  generally-used  boundary  condition,   but  for  small  values 
of  vi   and  large  \ ,    the  error  in  (2-99)   can  be  appreciable.     From 
(2-71)  we  obtain? 

»L,  *  "n.  +  ^n   *  "„/'-")  (2-100) 

where  the  subscript  /^refers  to  evaluation  at  a  distance  A^from  the 
wall. 

The  importance  of  equation  (2-71)  and  of  the  transformation  of 
the  continuity  equations  to  the  forms  (2-61)  and  (2-62)  is  now 
readily  seen.   The  normal  component  of  the  vector  u  defined  by 
equation  (2-60)  must  have  a  constant  magnitude,  --^j-  ,  along  the  sheath 
edge.  If  there  is  more  than  one  wall,  and  the  walls  involved  are 
orthogonal,  the  density  gradients  established  to  satisfy  this  condition 
at  each  sheath  edge  will  not  interfere  with  each  other.  Parallel 
surfaces,  on  the  other  hand,  create  an  eigenvalue  problem  because  u 
must  satisfy  the  boundary  condition  at  one  sheath  edge,  must  change 
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as  prescribed  by  the  continuity  equation  (2-6l)  as  we  move  across 
the  plasma,  and  must  acain  satisfy  uniquely  the  boundary  condition 
at  the  other  surface..  If  the  walls  are  neither  orthogonal  nor  parallel, 
the  density  gradients  established  to  meet  the  boundary  condition  for 
one  surface  will  have  along  the  other  surface  normal  components  which 

are  not  constant,  and  vice  versa.  This  will  result  in  a  violation 
of  equation  (2=71)  ■>  Consequently,  for  non-orthogonal  geometries  we 
are  unable  to  predict  the  steady  state  distributions  with  the  equations 
developed  in  this  thesis  except  by  an  approximate  method  to  be 
discussed  in  a  later  section. 
F.  Heat  Losses 

The  only  source  of  heat  loss  considered  in  this  thesis  is  that  which 
occurs  at  an  insulated  surface  as  a  consequence  of  the  particles 
impinging  on  the  surface  and  giving  up  their  energy.  This  energy  can 
come  from  several  souraess  (1)  the  recombination  energy,  eVj,  from 
each  electron-ion  recombination,  (2)  the  kinetic  energy  of  the  particles^ 
and  (3)  excitation  energy  of  excited  ions  which  return  to  their  stable 
state  upon  colliding  with  a  wall.  All  of  these  are  proportional 
to  the  particle  current  density  to  the  surfaces.  As  a  simplification 
for  illustrative  purposes  in  later  calculations  we  shall  limit  ourselves 
to  plasmas  whose  temperatures  are  sufficiently  low  (below  about  1CP 
degrees)  that  only  the  recombination  and  electron  kinetic  energies 
are  appreciable,.  The  heat  liberated  at  a  surface  is  then 

Eooall  -  Q  (eVl+   -i)j  (2-101) 

where  the  subscript  "n"  means  normal  to  the  surface  and  E.^^  is 

expressed  in  electron  volts. 
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Ill 

SOLUTIONS  OF  THE  CONTINUITY  EQUATION 

A.  General  Solutions 

This  section  discusses  a  procedure  to  obtain  particle  distri- 
bution functions  in  several  coordinate  systems,  and  presents  some 
of  the  general  features  of  the  solutions.  The  following  sections 
develop  many  solutions  for  specific  geometries. 

1.  Exact  Solutions 

a.  General  Developments 

With  no  magnetic  field,  we  have  exact  solutions  in  the 
eleven  orthogonal  coordinate  systems  in  which  the  wave  equation  is 
separable  when  the  boundary  surfaces  coincide  with  coordinate  planes. 
The  procedure  we  shall  describe  applies  to  the  non-steady  state 
plasma  for  vhich  we  repeat  the  general  continuity  equation  (2-£5) 


which  becomes  a  wave  equation  for  the  steady-state 


VV.         jL 


>,        %  s  (3-la> 

for  a  decaying  plasma  with  no  ionization  it  is 


RVV  -  §^  o-ib) 

with  a  steady°state  counterpart  (Laplace's  equation) 


7V  =  O.  O-lc) 

In  the  presence  of  a  magnetic  field,  we  further  require  that  the 

26 


coordinate  system  be  one  in  vhich  the  continuity  equation  may  be  ex- 
panded as 

Dm(1V?,V-f  T)r*Vj*.+  V7l.  --§£-  (3-2) 

(see  equation  (2-£7), 

Equations  (3=1)  and  (3-2)  are  separable  in  the  geometrical  coord- 
inates and  time,  and  the  general  solution  for  a  generalized  coordinate 
systems  u,  v,  w,  may  be  expressed  as 

n(u,v,w,t)  -  N(u,v,w,)T(t)„  (3-3) 

Substitution  of  equation  (3-3)  in  equation  (3-2),  and  division  by  NT 
yields 

D^+a^/+y=r=A         (M) 

where  -Kj.  is  a  separation  constant,  The  time  solution  is 


T*C,e-Kt  0-5) 


and  then 


In  the  absence  of  a  magnetic  field,  equation  (3-6)   reduces  to 

Vr-   -  -  k  ,  (3=7) 


in  which 


\ettiJ 


'*    y 


and,  for  the  steady  state  (k^  *  0) 

k1 "-   IX 
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Our  approach  to  the  solution  of  equations  (3=6)  and  (3=7)  is  to 
employ  the  technique  of  separation  of  variables,,  Thus  we  assume 
that  the  solution  is  of  the  form 

N(u,v,w)  -  f1(u)f2(v)f3(w))  (3=8) 

and,  by  substitution  into  equation  (3-6)  or  (3-7)  (as  appropriate), 

we  obtain  three  differential  equations,  one  in  terms  of  each  of  the 
coordinate  variables,  which  the  functions  f^,  f^s   and  ^must  separately 
satisfy,  Magnus  and  Oberhettinger  £9]  have  summarized  the  results  for 
eight  coordinate  systems,,  Because  the  separated  equations  are  of 
the  second  order,  there  will  be  two  linearly  independent  sets  of 
eigenf unctions  for  each  coordinate,  F'*'(kj,u)  and  F^2'(k^  ,u), 
G^!%jv)  and  G(2'(l%,v)$  and  H^  '( kn,w)  and  H^  '(kn,w),  which  satisfy 
the  differential  equation  in  the  corresponding  variable.  We  shall  call 
the  functions  with  the  superscript  (1)  functions  of  the  first  kind 
and  the  others  functions  of  the  second  kind.   The  most  general  solu= 
tion  to  these  equations  will  then  be  a  sum  of  the  products  of  all 
such  linearly  independent  solutions  of  the  form 

The  complete  solution  for  the  particle  density  distribution  is 


then 


n(«,Wt)  =  f~  ty„J«M»)Tfr)  (3_10) 
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The  summation  is  over  all  possible  values  of  jts  m,  and  n.  k^,  k^ 
and  kn  are  eigenvalues  or  spearation  constants  (depending  on  the 
coordinate  system  chosen)  corresponding  to  X,  m5  and  n  respectively*. 
VJe  can  sum  over  only  those  eigenf unctions  which  belong  to  a  degenerated 
eigenvalue  (as  determined  by  the  boundary  conditions).  It  should  be 
noted  that  the  time~dependent  solutions  to  be  developed  herein  are 
for  those  problems  for  which  the  boundary  condition  (2-71)  or  (2-99) 
is  still  applicable,,  A  plasma  decaying  from  a  steady-state  distri- 
bution upon  cessation  of  ionization  is  considered  to  be  such  a 
plasma. 

b.  Cartesian  Coordinates 

An  example  of  the  above  method  of  solution  is  the  solution 
in  cartesian  coordinates  in  which  the  functions  of  the  first  and 
second  kinds  may  all  be  cosines  and  sines  respectively^  or  in  which 
one  or  two  pairs  of  the  functions  may  be  exponentials.  For  examples 


y 

(3-U) 


By  substitution  of  this  solution  in  equation  (3-7)  >  we  see  that 
k  of  equations  (3~7a)  and  (3»7b)  is  given  by 

and  consequently  that  the  time  solution  (equation  (3=5))  is 


**% 
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Obviously,  when 

the  tine  dependence  vanishes,  and  the  plasma  is  "stationary,," 
Another  possible  solution,   involving  exponentials,   is 

Nty •>  -£, fc/'^e ^^fe^H^^^^  ^7**^4 (3-13) 
where  k2  is  given  by  k^  =    Z    (£>2-^/-^)^  (3-13a) 

c.     Cylindrical  Coordinates 

There  is  an  independent  eigenvalue,  k«,  for  each  set  of 
functions  only  in  cartesian  coordinates.     For  any  other  coordinate 
system  u,v,w,  defined  by  a  transformation  equation  of  the  form 

x  «  x(u,v,w) 

y  s  y(u,v,w)  (3-lU) 

z  ■  z(u,v,w). 
those  eigenf unctions  of  the  variables  u,v,w  which  appear  together  in 
any  of  equations  (3-lh)  will  be  connected  through  a  separation  variable 
which  is  required  to  reduce  equation  (3-li)  or  (3-7)  to  separate 
equations  in  each  of  the  coordinates,  Each  such  separation  variable 
required  for  the  above  purpose  reduces  the  number  of  eigenvalues 
available  for  the  satisfaction  of  the  boundary  conditions.  For  ex= 
ample,  for  cylindrical  coordinates,  where 

x  ■  recoil 

y  ■  r&o*1? 

z  -  z,  (3-15) 

the  eigenf  unctions  in   r  and  l(>  are  related  by  the  separation  constant, 
m,  as  indicated  in  the  we  Unknown  solution 
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(3-16) 


vhere 


and 


"4*  +  C7, 


feH=<^ 


k2--   16r/*C2  (3-l6b) 

An  alternative  Z(z)  solution  is 

^.^e1'1^/;^  i„yo        (3-i6c) 

vith 


( 3-l6d) 


JpjCk  r)  and  \(k  r)  are  Bessel  functions  of  order  m  of  the  first  and 
second  kinds  respectively  (Bessel  and  Neumann  functions). 

c.  Spherical^  Prolate  and  Oblate  Ellipsoidal  Coordinates 

In  the  other  coordinate  systems  to  be  discussed,  the  general 
solutions  of  (3-7)  ares 
Spherical s 
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(3-17) 


vrhere 


and 


(3-17b) 


and  in  which  pf5(cos©)  and  Qn(eosQ)  are  Legendre  Polynomials  of  the  first 
and  second  kinds  respectively  and  jn(k  r)  and  nn(k  r)  are  spherical 
Bessel  functions  of  order  n  of  the  first  and  second  kinds  respectively* 
The  latter  are  defined  by  the  equations 

If 


Prolate  Ellipsoidal  (See  Appendix  I  for  diagram  and  transformation 
equations!  notation  employed  is  that  of  Flammer  £p3  with  slight  modi- 
fication) % 

au  (f  if) * &**$£  (cy  t  cilm„  s2  Ml  * 

where 


and  d  is  the  interfocal  distance,  S^  (c,  9)  and  sjjp^ ,  ^)  are  spheroidal 
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angle    functions  of  the  first  and  second  kinds  and  R^n '(cj    \)   and 
r(2)(  a      2'  )  are  spheroidal  radial  functions.  See  Appendix  II  for 

tables  and  graphs  of  these  functions. 

Oblate  Ellipsoidal  functions  may  be  obtained  from  the  prolate  func- 
tions by  the  substitution  of  =ic,  for  c*  where  appearing  in  equation 
(3-18).  Equation  (3-l8a)  is  unchanged  by  the  transformation. 

Similar  equations  are  obtained  in  parabolic  cylinder  coordinates 
(where  the  functions  include  Weber-Hermite  functions),  parabolic 
coordinates  (Vhittaker  functions),  and  elliptic  cylinder  coordinates 
(Mathieu  functions).  These  are  not  developed  further  in  this  thesis 
because  their  geometries  are  not  useful  in  the  problems  with  vhich  we 
shall  deal. 

d.  Solutions  in  the  Presence  of  a  Magnetic  Field 

In  the  presence  of  a  magnetic  field,  Dm  is  a  tensor,  and 
equation  (3-6)  must  be  employed.  In  order  to  make  a  separation  of 
the  form  (3-6)>  we  are  limited  to  cartesian  coordinates  or  cylindrical 
coordinates  for  the  problems  which  we  shall  treat.  Parabolic  cylinder 
and  elliptic  cylinder  coordinates  are  also  suitable,  but  we  shall  not 
need  them  for  the  geometries  with  which  we  shall  be  concerned.  Ftrther, 
if  we  consider  the  assumptions  under  which  T)m   was  developed,  we  may 
immediately  limit  the  treatment  to  the  case  where  m  and  C  are  zero, 
since  no  gradients  in  the  $  direction  can  exist. 

In  the  coordinate  systems  selected,  equation  (3-6)  can 
be  reduced  by  separation  of  variables,  by  assuming 
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and  thus  obtaining  two  equations 


^X  „  _  ¥ 


Ml     K 


and 


(3-20a) 


V'^"  =  -J&L   -.  .I'''  (3-20b) 


where 


V«l+l^A*"k,-|!*^-  (3"20c> 

The  solution  to  the  above  equations  is  the  sane  as  in  the  absence  of  a 
field,  except  that  equations  (3-12)  and  (3-l6d)  become  (for  a  field  in 
the  z-direction) s 

and 


c  f         i  ?      T>      U  *-"\  (3-l6d«) 

e.  Further  Simplification  of  Previous  Solutions 

The  equations  developed  so  far  can  be  simplified  considerably 
by  considering  some  of  the  restrictions  upon  the  solutions,,  The  require- 
ment  of  a  single-valuedness  limits,/,  m,  and  n  to  be  integers,  and  it 
is  possible  to  further  limit  them  to  positive  values  without  loss  of 
generality.  Further,  if  the  plasma  possesses  rotational  symmetry 

(which  it  must  in  the  presence  of  a  magnetic  field  under  our  assump- 
tions), then  m  and  C5  are  both  zero  in   the  equations  for  I\I(u,v,w) 
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(vith  the  exception  of  rectangular  coordinates).  The  requirement 
that  the  density  be  non-negative  everywhere  within  a  given  region  vill 
further  limit  the  solutions.  Because  all  the  functions  discussed 
(except  the  exponential  function)  are  periodic,  this  latter  requirement 
usually  limits  the  steady-state  solutions  for  bounded  regions  to  the 
first  positive  segment  of  the  applicable  function,  and  values  of  k, , 
kpj,  and/or  kn  are  thus  determined  by  the  dimensions  of  the  plasma 
container.  From  the  boundary  condition  (2-71),  the  physical  reason  for 
this  effect  is  apparent.   u  ^-j-has  a  definite  value  at  the  edge  of 
the  sheath.  Equation  (3-7)  can  be  transformed  into  an  equation  in  u 
in  the  same  manner  as  in  the  derivation  of  equation  (2-62),  yielding 

V-U   +  u"*  e-k*  (3-21) 

It  is  the  value  of  kc   in  this  equation  that  determines  the  "trajec- 
tory" of  the  vector  u  across  the  plasma  to  the  sheath  en  the  "opposite" 
side  where  u  must  again  satisfy  the  boundary  condition  (2-71).  Because 
equation  (3-22)  is  a  first  order  equation  in  u,  the  problem  is  completely 
stated  (except  for  a  condition  to  determine  an  arbitrary  normalization 
constant),  and  only  one  value  of  k  is  possible. 

Further.,  and  obvious 5  simplification  can  be  accomplished  for  many 
of  the  solutions  stated  herein  if  they  must  describe  a  plasma  in  a 

region  in  thich  some  of  the  eigenfunctions  are  singular.  In  parti- 
cular, the  following  functions  have  singularities  for  the  values  indicated; 
%CO),f|n(0),  QS(+1),  4ftc,il),  S&2)  (-ic,!l). 
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2.  Approximate  Solutions 

There  are  a  limited  number  of  geometries  for  vhich  a  stable 
ambipolar  solution  may  be  obtained  by  the  separation  procedure  de- 
scribed herein,,  It  is,  however,  possible  to  obtain  an  approximate 
solution  for  the  particle  density  distribution  if  the  geometry  is 
not  too  radically  different  from  a  geometry  for  which  an  exact 
solution  iscbtainableo 

In  the  discussion  of  the  boundary  conditions  (Section  II-E-3)> 
it  was  shown  that  the  conditions  at  the  sheath  edge  control  the  den- 
sity distribution  in  the  vicinity  of  a  surface.  If  there  are  two  or 
more  surfaces  bounding  the  plasma  in  such  a  manner  that  it  is  not 
possible  to  obtain  an  exact  solution  for  the  particle  density  dis- 
tribution, then  we  divide  the  plasma  into  "zones  of  influence" 
about  each  of  these  surfaces .  We  then  calculate  the  density  dis- 
tribution for  each  zone  as  determined  by  the  boundary  conditions  at 
the  insulator  surfaces  belonging  to  that  zone.  In  the  region  near 
such  surfaces  we  expect  our  solution  to  be  very  nearly  correct.  There 
will  be  a  transition  region  between  the  zones  of  influence  in  which 
we  expect  the  greatest  difference  between  the  model  and  the  actual 
plasma.  In  this  region  there  is  either  a  departure  from  the  ambi- 
polar diffusion  or  a  time  dependence  (non-stationary  plasma),  or 
both.  However,  we  expect  that  processes  in  this  transition 
region  do  not  affect  the  total  solutions  significantly  in  many 
problems.  A  measure  of  the  importance 
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of  these  processes  Is  the  degree  to  vihich  the  solutions  for  the 
different  zones  conflict  in  the  transition  region.  When  the  dis- 
continuity across  the  transition  region  is  small,  we  expect  that  our 
solution  correctly  describes  the  plasma  in  the  immediate  vicinity 
of  the  insulator  walls  and  describes  it  approximately  (within  an 
order  of  magnitude)  even  in  the  transition  region. 

In  order  to  discuss  this  approximation  more  quantitatively,  let 
us  consider  two  walls,  W^  and  Wp>  bounding  a  partly- ionized  gas,  and 
assume  that  the  configurations  are  such  that  the  model  described  in 
this  thesis  can  be  applied  to  each  wall  in  the  absence  of  the  other. 
Thus,  two  distribution  functions,  n,(q.)  and  ^(qj),  are  obtained 
(with  qj  being  the  coordinates).  If,  through  adjustment  of  the 
spearation  variables  of  equations  (3-6)  and  (3-7)  and  of  the  normali- 
zation constants,  it  is  possible  to  make 

T?,^*  (3-2lia) 

and 

everywhere,  then  the  solution  is  the  exact  solution.  If  this  cannot 
be  done,  then  we  must  find  a  surface,  S,  which  separates  the  zones 
of  influence  of  the  walls  W^  and  ¥2.  The  separate  solutions  for 
these  zones,  joined  at  S,  form  the  zero~order  approximation  of  our 
solution. 

There  are  three  general  procedures  by  which  we  might  develop  a 
criterion  for  determining  the  position  of  the  surface  S.  We  shall 
briefly  discuss  all  three  in  the  following  paragraphs,  but  shall  not 
attempt  a  rigorous  proof. 
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One  procedure  is  to  employ  a  perturbation  theory  in  which  we 
define 

71±--riG-t   ^7?t'-»-.  ..  (3-25) 

and 

E=E0+MMV.,.  (3-26) 

where  ){    is  a  parameter  defined  by  the  equation 


??O0V'E  -  ^fa.-*-) 


(3-27) 


in  which  n_0   is  defined  as   the  maximum  density  in  the  plasma.     From 
equation  (3-27)  and  Poisson*s  equation  it   is  apparent  that 

* =  it  (3~28) 

For  our  zero  approximation,  we  impose  the  conditions! 

and 

(L    -P-0-P6  »-30) 

in  each  zone  of  influence. 

In  order  to  determine  a  criterion  for  positioning  the  surface 
between  the  zones  of  influence,  S,   (i.e.,   in  order  to  find  an  adapted 
zero-order  approximation),  we  must  look  at  the  first  order  terms. 
Neglecting  quadratic  terms  in  )f\  ,  we  obtain 


PJ  =  /a,(eQ  V  +  no  E')  -  D+VtiJ  (3-3U) 

pj  =      -/l/fo*-'  <"  *  *'  )  -   D"^7?:  (3"31b) 
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If  we  further  assume  that  either  there  is  no  time  dependence  or  that 
the  time  dependence  can  be  accounted  for  by  replacing  all  quantities 
in  the  following  equations  by  their  time-averages,   we  obtain  the 
following  from  the  continuity  equation: 

WT  -tm'  t  $(?-& *  •fiV'JW)  +  ^(r-r%)(ai  - f%  )  (3-38.) 

and 

where  the  subscript  "S"  means  evaluation  at  the  surface  S,   ^(r-r<.N) 
is  the  Dirac  delta  function,  and  n  is  the  unit  normal  to  the  surface 
S.  The  last  two  terms  in  the  above  equations  represent  particles 
which  would  have  to  be  created  or  absorbed  at  the  imaginary  surface 
S  in  order  to  supply  the  currents  indicated.  The  first  order  density 
terms  can  be  either  positive  or  negative.  It  is  apparent  from  (3-32) 
that  at  the  surface  S, 

jifo,, -?W)  •-  -(K  -K>*-&.  -^)      (3-33.) 

If  the  terms  of  equation  (3-33)  are  identically  zero,  then  the 
first  order  solution  can  be  the1  trivial  solution 

71±  =0  (3-31*) 

Any  deviation  of  the  first  order  solution  from  equation  C 3—3U)  is 
obviously  a  function  of  the  magnitude  of  the  first  term  in  equation 
(3-33).  The  surface  S  is  thus  the  surface  at  which  this  term  is  a 
minimum.  If  there  were  a  surface  along  which 

(Dm*0H  -(DM'"***).*  (3"35) 
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everywhere,   subject  to  the  condition  that  the  density  solutions  were 
also  equal  at  the  surface,   then  the       iJf^o\,  terms  would 

vanish.     This   is  not  possible  generally,  however,  and  a  least-mean- 
square  procedure  should  be  applied.     Because  this   is  not  feasible 
for  many  problems,   the   following  criterion  for  determining  the 
position  of  S,  applicable  at  a  given  point  in  which  we  are   interested, 
is  employed  to  obtain  our  adjusted  zero~order  approximations 
|  Dm    Vn0|o5or,|    =  I0™'V*02  «o|kl 

where     Y  \   is   the  angle  between    Vf\ot  and  the  surface  normal  and     Xx 

_^ 
between  V/iOJr  and  the  surface  normal.       Further,    if 


?7«  -  Cfii  ^ J 


(3-37a) 


and 


then  at  the  surface  S 

Solving  for  C2,  substituting  the  result  in  equation  (3-36),  and 
dividing  both  sides  by  C^f^(qj)s  yields 


(3-37b) 


(3-38) 


D 


fY\ 


ccz> 


r. 


(3-39) 


Because  only  order  of  magnitude  results  will  be  obtainable  from 
this  approximation,   it  is  reasonable  to  assume  as  a  first  approxima- 
tion that 

OStf,    *   cod  Ik  (3-UO) 

ho 


By  employing  equations   (3-37),  ve  obtain  the  following  conditions 
defining  the  surface  of  separation 


The  second  procedure  to  derive  equation   ( 3-Ul)    is  to  consider 
the  relative   influence  of  each  boundary  at  all  points  within  the 
plasma,  and  define  the  surface  S  as  that  surface  at  which  the   in- 
fluences of  the  two  boundaries  are  equal.     Some  measure  of  this 
influence  must  first  be  found.     If  we  have   in  the  neighborhood  of 
a  boundary  a  solution  for  the  particle  density  which  has  a  maximum, 
then  certainly  beyond  that  maximum  the  plasma   is  not    influenced 
appreciably  by  that  boundary.     For  the  region  between  the  surface 
and  the  maximum,   the  measure  of  a  boundary's   influence  could  be 
the  current  created  in  response  to  its  presence.     This  procedure 
leads  to  equation  (3-^1)   with  an  added  restriction  that   if  there 
are  two  points  at  which  this  equation  is  satisfied,  we  must  select 
the  one  within  the  zones  defined  by  the  walls  and  the  first  maxima 
of  their  respective  solutions  whenever  possible. 

The   third  procedure  to  determine   the  position  of  the  surface 
S  would  be  to  calculate  the  probability   that  a  particle  at  a  certain 
position  would  go  to  the  one  wall  or  the  other.     The  surface  S 
dividing   the  two  zones  would  be  the  surface  along  which  particles 
would  have  equal  probabilities  of  reaching  all  the  walls   in  question. 
This  might  be  calculated  by  stochastic  methods  or  from  the  diffusion 
equations.     We  did  not  attempt  the  former.     The   latter  method   is 
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to  find  the  surface  from  which  the  mathematical  solution  indicates 
that  as  many  of  the  particles  produced  within  a  differential  volume 
straddling  the  surface  would  go  to  the  one  zone  as  to  the  other. 
This  procedure  also  leads  to  equation  (3-lU)  as  the  criterion. 
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B.     Specific  Solutions  of  the  Continuity  Equations. 
1.     General  Conditions 

In  order  to  obtain  an  estimate  of  the  magnitude  of  certain 
quantities   in  the  equations  to  follow,   we  shall  use  the  following 
sample  data,  some  of  which  is  taken  from  Brown  [1/]  % 
gass       Helium 
p     e     1  mm  Hg« 
T     -     700°K 

T     =     50,000°K  *s>U.3  eV. 
ft+    ^s   10^  cm^/volt-sec    as  1  m2/volt-sec 
/U.  fin  7»6  X  ic£     cm2/volt-sec  =?  76  m2/volt-sec  ' 

D+  $&     ^ffa        **  «°6  n  /sec 

D.  -s      SE  ^ _        se»         3U0  n2/sec 
Ds   c*  Da  -     q/u-+P-/v  =  U.U  m2/sec 

,,  ...         -^       lull  volts 
eVi    ^»    2U.58  eV 

eVi  +   1?    ~       2?  eV 
E     ■  300  volts/m 

n+Q    ^     n„0     -  r^     »  1012  cm"3  ■  10l8m~3 

where  "r^1  is  the  density  in  the  center  of  the  plasma  or  in  an  un- 
disturbed, constant-density  region9 

For  the  calculations  of  this  section  we  assume  that  the  mean  free 
paths  of  the  ions  and  electrons  are  both  sufficiently  smaller  than  any 
boundary  dimensions   involved,  so  that  boundary  condition  (2-99), 

Km  -'  °- 
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is  applicable.  The  more  general  boundary  condition  (2-71)  will  be 
employed  in  the  discussion  of  cylindrical  probes  in  the  next  chapter,, 
and  the  method  of  employment  of  that  boundary  condition  will  be  dis- 
cussed then. 

A  quantity  often  quoted  and  employed  in  the  literature  is  the 
diffusion  length, J\.s   which  is  defined  by  the  equation 

We  shall  compute  this  quantity  for  the  various  geometries  in  this 
section. 

Finally,  we  shall  calculate  the  "zone  of  influence"  of  many  of 
the  surfaces  to  be  discussed.  This  we  define  as  that  region  in  which 
the  particle  density  is  less  than  £(#  or  less  than  90$  (as  indicated) 
of  the  density  at  the  center  or  undisturbed  region  of  the  plasma. 

2.  Plane  and  Rectangular  Configurations 
a.  Rectangular  cavity 

Because  it  is  the  most  complex  geometry  we  can  treat  in 
rectangular  coordinates,  we  begin  with  a  discussion  of  a  rectangular 
cavity  of  dimensions  Xq,  y0,  and  z0.  The  limitations  of  the  develop- 
ment of  the  diffusion  tensor  prevent  the  application  of  the  follow- 
ing solutions  to  magnetic  fields,  and  thus  the  scalar  coefficient  Ds 
will  be  used.  Of  the  possible  variations  of  the  solution  (such  as 
equations  (3-11)  and  (3-13)),  we  can  immediately  rule  out  those  con- 
taining the  exponentials  because  of  the  requirement  of  zero  density 

at  the  walls.  The  boundary  conditions  also  limit  the  solution  to 
the  first  positive  segments  of  the  sine  functions  (and  exclude  the 
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cosine  functions)  if  the  origin  is  taken  at  one  corner  of  the  cavity. 
If  the  origin  vere  selected  at  the  center  of  the  cavity,  the  cosine 
functions  would  be  used  and  the  sine  functions  excluded.  With  the 
origin  at  one  corner,  and  the  sides  along  the  positive  coordinate  axes, 
we  obtain  the  following  solution  for  the  steady  state  (k^  *  0)  plasma; 

Ttfty*)   -  yio^i\^)^4\{)y^.(l,i)  (3-52) 

with  the  diffusion  length  given  by 

3P         *T  "  Ds  "  (*~*>     ^'      ta*'  (3"52a) 

The  three  components  of  the  particle  current  density  at  any  point 
are  easily  obtained,  using  equations  (2-li2) 

P„  --  -rfc  Ds  I  to*(l*)  *Xfa)  W4*)       (3-53) 

with  similar  equations  for  the  other  components.  From  these  equations 
it  is  seen  that  the  particle  current  density  is  zero  at  the  center  of 
the  cavity,  is  directed  outward  from  the  center  to  the  walls,  and  is 
a  maximum  at  the  walls  where,  for  the  data  of  (3-5>0)  and  a  £  centi- 
meter wall-to-wall  separation,  the  magnitude  of  the  current  density 
at  the  center  of  a  wall  would  be  approximately  2.7  X  10   particles/cm^ 
-sec.  It  is  of  interest  to  note  that  the  steady -state  particle  current 
densities  at  the  walls  decrease  in  inverse  proportion  to  a  linear  in- 
crease in  cavity  dimensions 0  This  decrease  is  a  consequence  of  the 
decrease  in  the  gradient  of  the  current  density.  From  equation  (3-!?2a) 
one  can  see  that  the  rate  of  ionization,  ))     ,  necessary  to  replace  the 
particle  losses  to  the  walls  is  inversely  proportional  to  the  square 
of  the  dimensions  (and,  incidentally,  is  primarily  determined  by  the 
smallest  dimension  of 

16 


the  cavity  as  a  consequence),,  The  explanation  for  the  inverse 
dependence  of  the  rate  of  ionization  per  unit  volume  is  easily  seen 
by  equating  the  number  of  particles  destroyed  at  the  vails  to  the 
total  number  created  within  the  cavity.  The  latter  is  proportional 
to  the  volume  and  to  )}    0  The  number  of  particles  carried,  to  the  walls 
however,  is  proportional  to  the  surface  of  the  walls  and  to  the  current 
density.  The  latter  is  inversely  proportional  to  the  linear  dimen- 
sion, as  previously  mentioned,  and  consequently  the  number  of  parti- 
cles  drawn  off  is  proportional  only  to  the  first  power  of  the  linear 
dimension  of  the  cavity.  Consequently,  the  inverse  dependent  of 
on  the  linear  dimensions  is  expected.  The  fact  that  the  number  of 
particles  destroyed  at  the  walls  is  directly  proportional  to  the 
linear  dimensions  of  the  cavity  (and  therefore,  so  are  the  energy 
losses)  is  also  worth  noting,, 

The  space  charge  field,  according  to  equation  (2-38)  is 

or,  again  using  the  data  from  (3-£0),  is  approximately  =270  volts/m 
(directed  away  from  the  center  of  the  cavity)  at  the  point  (.2£xQ, 
•5y0?  °£z0).  The  cavity  dimensions  are  assumed  to  be  £cn.  Again, 
an  inverse  relationship  to  the  wall-to-wall  separation  exists.  The 
space  charge  potential  (equation  (2»Iil))  is 

or  approximately  =1.53  volts  at  the  point  (.25x  ,  .5y  ,  «5Z0)  for 
a£  cm.  cavity.  The  potential  is  with  respect  to  the  center  of  the 
cavity.  The  potential  predicted  by  equation  (3-5h)   becomes  singular 
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at  the  vails,  and  this  feature  of  the  solution  indicates  a  limitation 
of  these  equations 2  they  are  not  applicable  beyond  the  sheath  edge. 

Using  the  current  equations  (3-53)  we  can  compute  the  energy 
lost  at  the  walls.  From  equation  (2-101)  ve  obtain 

where  the  factor  2  takes  into  account  the  surfaces  adjoining  the  corner 
opposite  the  origin.  This  equation  reduces  to 


(3-5*0 


E^i ? fev  **+*£] 

271oDs(e V  *£)  f^^ *  *****  +  ^V] 
The  heat  liberated  per  unit  area  is  obtained  by  dividing  each  integral 
in  (3-55)  by  the  appropriate  area,  obtaining  for  the  three  sides  of 
the  cavity: 

/TV 

sfctv*oT  TXc  "  >  "^  )  T*o         ^(3-55b) 

Again  using  5  cm.  dimensions  and  the  data  from  (3-50),  we  see  that 
EVTall  is  approximately  1^.86  X  1019  ev/sec^7.8  joules/sec.     From 
equation  (3-55a)    it  is  seen  again  that  the  total  heat  lost  per  secod 
is  proportional  to  the  linear  dimensions,  and  that  doubling  the 
length  of  all  the  sides,  for  instance,  will  double  the  rate  of  heat 
loss. 

One  very  important  fact  about  plasma  characteristics  may  be  seen 
from  the  preceeding  equations,  and  in  particular  from  equation  (3-52a). 
That  is  that  the  effect  of  a  small  dimension  on  the  characteristics 
is  very  much  greater  than  that  of  large  dimensions,   so  that  a  long, 
narrow  plasma  has,   in  the  middle  region  of  the  longer  dimension, 
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characteristics  that  are  almost  completely  dependent  on  the  shorter 
dimension  and  are  relatively  independent  of  the  longer  dimension.  Thus, 
a  plasma  whose  z-dimension  is  very  much  greater  than  its  dimensions 
in  other  directions  can  often  be  treated  quite  properly  as  "infinite" 
in  the  z~d£reet£cn,  with  corresponding  simplifications  of  the  equa- 
tions describing  the  plasma. 

The  zones  of  influence  of  the  rectangular  surfaces  are  indicated 
in  Figure  3-1 • 

If  ionization  ceases,  ))  «  0,  and  the  plasma  decays  exponentially 
in  accordance  with  equation  (3=5),  with  kt  given  by  (3-7a),  (3-12), 
and  (3-£2a)  as 

K-(^(l)M£t  (3-56) 

Thus  the  plasma  will  decay  to  l/e  of  its  initial  value  in  time. 

which,  for  a  cavity  of  £  cm  dimensions,  is  about  20  Usee  This 
decay  time  is  inversely  proportional  to  the  square  of  the  linear 
dimensions. 

If  the  iritial  distribution  is  not  the  steady-state  distrib 
tior,  but  instead  some  function  f(x,y,z),  then  the  coefficients  of 
equation  (3=11)  must  be  obtained  in  the  normal  manner  for  a  Fourier 
series,  and  a  series  of  the  form  (3-10)  must  be  developed.  For  example, 
if  the  initial  distribution  is  a  constant  density  n  throughout 
(except  discontinuous ly  zero  at  the  surfaces),  the  particle  distribution 
function  will  be 

1        ''""""  1,8  (3-58) 


1x9 


with 


h 


(3-58a) 


It  is  to  be  noted  that  the  higher  order  terms  decay  at  a  faster  rate, 
so  that  the  plasma  approaches  the  steady-state  distribution  as  it  decays. 
Figure  3-2  shows  the  decay  of  such  a  plasma  at  three  different  times. 
Other  features  of  this  solution  are  that  the  rate  of  decay  is  in- 
versely  proportional  to  the  square  of  the  linear  dimensions  so  that 
a  plasma  in  a  large  cavity  decays  more  slowly  than  one  in  a  small  one. 
Further,  we  again  see  that  the  shortest  dimensions  control  the  decay 
rate,  so  that  dimensions  which  are  very  much  larger  than  the  smallest 
dimension  may  again  be  treated  as  infinite.  For  instance,  for  a 
length  which  is  ten  times  that  of  the. others,  the  error  in  ignoring 
it  in  equation  (3-57)  is  less  than  one  percent. 

b.  Infinite  tubes  rectangular  cross-section 

By  allowing  zQ   to  go  to  infinity  and  setting    %  s  y» 
in  the  equations  of  the  previous  section,  we  obtain  the  equations  for 
an  infinitely  long  rectangular  geometry.  Thus  we  have  for  the  steady 
state 

with 

The  current  density?  potential  and  other  equations  are  similarly  ob- 
tained and  will  not  be  stated  here. 
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c.  Infinite  Plane» 

By  again  allowing  one  dimension  to  become  infinite,  for 
example,  y0,  and  by  replacing  y  by  y0/2,  we  obtain  the  equations  for 
the  simplest  geometry,  two  parallel  planes  separated  by  a  distance  xQ, 
or  a  single  plane  in  contact  with  an  infinite  plasma  in  which 

x.  =  *^,  o-60* 

where  we  know  the  rate  of  ionization  from  other  considerations,  or 
where  x0/2  is  the  distanceto  the  undisturbed  region  of  the  plasma. 
The  restriction  that  there  be  no  magnetic  field  can  now  be  dropped  if 
the  only  density  gradient  is  along  the  coordinate  lines  perpendicular 
to  the  planes  and  if  the  field  is  parallel  to  these  planes.  If  the 
field  is  perpendicular  to  the  plane  (so  that  it  does  not  affect  diffusion 
perpendicular  to  the  plane) ,  or  if  there  is  no  field,  then  Dmi  in 
equation  (3-60)  and  in  the  equations  to  follow  can  be  replaced  by  its 
no-field  limit,  Ds„ 

The  particle  density  distribution  for  the  steady  state  is 

7t'7l**J*.($ix),  (3-61) 

with 

■t'V'A'   ftV"'  (3-61*) 

The  other  equations  are  derived  from  those  describing  the  rectangular 
cavity  by  replacing  y  and  z  by  y0/2  and  z0/2  respectively,  by  letting 
y0  and  z0  become  infinite,  and  by  replacing  Ds  by  Dm,  . 

3.  Cylindrical  Configurations 
a.  Finite  Cylinder,  interior 

For  a  cylinder  of  radius  R  and  length  zQ,  we  can  immediately 
simplify  Equation  (3=16)  for  the  steady  state  by  considering  the 
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functions  and  the  boundary  conditions.  The  Neumann  functions  are 
singular  at  the  origin*  and  so  their  coefficients,  Cofo*  must  De 
zero.  Further,  if  the  origin  is  taken  at  the  center  of  the  Ifase  of 
the  cylinder,  the  coefficients  of  the  cosine  functions,  Cyn,  must  be 
zero  also,  because  the  density  is  zero  at  z  ■  0  (using  boundary 
condition  (2-99)).  Then,  as  in  the  case  of  the  rectangular  configura- 
tions, the  requirements  that  n  always  be  positive  (or  zero)  and  the 
boundary  condition  at  z  ^  z0  requires  us  to  set  all  the  sine 
coefficients  to  zero  except  C^*,  and  kg  becomes  ^   •  Because  ve 
have  discussed  no  forces  vhich  would  tend  to  establish  an  assymetric 
angular  distribution,  and  indeed  in  the  case  of  the  derivation  of 
the  diffusion  tensor  ir  a  magnetic  field  have  specified  that  no 
such  assymetry  exist,  we  set  m  «  0.  Thus  we  are  limited  to  the 
zero  order  Bessel  function,  and  because  of  the  boundary  conditions 
must  select  the  first  positive  segment  of  that  function.  Thus  we  ob- 
tain the  well-known  distribution  function 

^^ni.(^WiA  (3"62) 


where 


J_    ±     4  J-   -/.LW+^Y"  (3-62a) 


and 


y=(^X+(i)V„  (>62b) 

The  magnetic  field,  if  any,  is  parallel  to  the  cylinder  axis. 
The  ion  current  density  components  are  obtained  by  using 
Equation  ( 2-1*2)  s 
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Thus,  at  the  center  of  the  cylinder  (r  =  0,  z  -   -|-  z0),  the  current 
density  is  zero,  and  within  the  remainder  of  the  cylinder  I  ,-   and  I  g 
are  both  directed  from  the  center  to  the  walls9   |V  *s  a  maximum 
at  the  ends  of  the  cylinder,  but  I  .  reaches.,  a  maximum  (see  Figure 
3-U)  when 

i,e*,  at  r  -  0,75>R»  This  maximum  occurs  at  a  radius  less  than  R 
because  the  additional  current  supplied  by  ionization  beyond  r  -  0,7!?R 
is  not  sufficient  to  counteract  the  effect  of  the  increasing  cylin- 
drical volume  into  which  the  current  is  flowing.  As  for  the  rectangu- 
lar cavity,  an  increase  in  cylinder  size  is  accompanied  by  a  decrease 
in  ionization  rate  and,  current  densities  in  the  steady  state. 

The  space  charge  fields  as  obtained  from  (2-35)  and  (2=36)  is: 

The  zones  of  influence  of  the  cylinder  surface  are  indicated 
in  Figure  3-3« 

One  of  the  characteristics  of  this  plasma  which  is  of  interest 
is  the  fact  that  the  shape  of  the  steady  state  configuration  of  ion 
density  is  dependent  only  upon  the  cylinder  dimensions  if  the  strength 
of  the  magnetic  field  is  less  than  that  at  which  the  "pinch"  effect 


£r  =  ,£j_  =  •- 


begins  to  appear.  However,  the  ion  current  density  in  the  radial 
direction  (perpendicular  to  the  field)  is  proportional  to  the  inverse 
square  of  the  magnetic  field  (at  high  B)  through  its  dependence  on 
Dm  .  . 

If  the  steady  state  plasma  begins  to  decay,  the  distribution 
is  given  by  the  equations 

*r,,t>-  v^(li)^r)em°"^mV^]\  (3-65) 

If,  however,  a  decaying  plasma  with  initial  distribution  F(r,l|,z) 
is  to  be  described,  it  is  necessary  to  develop  a  series  from  Equation 
(3-16)  in  the  same  manner  as  used  to  develop  a  Fourier  series,  obtain- 
ing equations  similar  to  (3-5>8) « 

b.  Infinite  Cylinder,  interior 

In  a  manner  similar  to  that  of  the  rectangular  problem 
it  is  possible  to  treat  a  very  long  cylinder  as  if  it  were  effectively 
of  infinite  length.  Employing  identical  procedures,  ve  obtain  for 
the  particle  density  distribution  function  of  a  "long"  cylinder  of 
radius  Rs 

H(r)  =  *o  X^'r),  (3-66) 

where 

I    /2,V\2-       n    u   /l^\2rs  (3-66a) 

The  particle  current  density  is  given  by 


Pr  -    +Z-j?  D„^Jt(^r)         H,  -  Pz  -  O,  (3-67) 


the  field  by 
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the  space  charge  potential  by 

v-£±£!^k^  (3.69) 

and  the  heat  liberated  at  the  walls  per  unit  cylinder  length  by 

/  ir\  ( 3-70) 

It  should  be  noted  that  the  latter  is  independent  of  the  cylinder 
dimensions,   and  using  the  data  of  (3-i>0),    is  approximately  1,32  X10  y 
g-2  eV/m-sec  -  2,11  B       joule/m~sec     (with  B  in  webers/m2). 

The  decay  problem  is  treated  in  the  same  manner  as  for  the  finite 
cylinder,  and  we  obtain  for  decay  from  the  steady  state 

x(n&  -  HM'-jbe't***7*?  (3-71) 

Figure  3-k   is  a  plot  of  the  radial  current  density  as  a  function 
of  radial  position  for  a  steady  plasma.  Figure  3-5  is  a  graph  of 
the  space  charge  potential  versus  radial  position  for  an  infinite 
cylinder.  The  fact  that  the  space  charge  field  and  potential  become 
infinite  at  the  wall  is  a  consequence  of  the  assumption  that  J\^n 
■  0.  As  in  the  rectangular  cases  these  equations  are  not  valid  within 
the  sheath. 

c*  Infinite  Cylinder  in  Infinite  Plasma 

Using  our  techniques,  there  is  no  exact  steady-state 
solution  to  the  continuity  equation  (3-1)  or  (3-2)  for  a  plasma  exterior 
to  a  cylinder  of  finite  length  because  we  cannot  combine  the  solutions 
(3-16)  in  a  manner  such  as  to  meet  the  boundary  conditions  at  the  ends 
of  the  cylinder.   However,  for  the  infinite  cylinder  of  radius  R 
it  is  possible  to  solve  for  the  distribution  function.  Again  we  limit 
ourselves  to  those  problems  in  which  no  assymetric  forces  are  acting 


(and  thus  to  axialiy  symmetric  plasmas),  but  retain  the  Neumann 
functions  because  the  origin  is  no  longer  within  the  limits  of  r.     VJe 
obtain  from  (3-16),  for  the  steady  state, 

*fr)  *  A  [U^)  *  ChJ0(}&Lr){         *i  r  *  *'  (3-72) 


vrhere  A  is  a  normalization  constant  which  we  shall  determine  later, 

R  is  the  cylinder  radius,  RB  is  that  radius  at  which  n(r)  is  a  maximum, 

and 

-X  (  \SL   r.  A 

(3-72a) 


c  =  - 

if  the  boundary  condition  (2-99;  ('Hjall  e  ^)  applies,  and  is  a  more 
complicated  function  if  equation  (2-71)  applies.  The  latter  function 
will  be  developed  in  the  section  on  probe  theory*  RJ'  the  radius  at 
which  equation  (3-72)  is  a  maximum,  is  found  by  setting  the  first 
derivative  to  zero;  i.e.,  it  is  the  solution  of 

The  normalization  constant,  A,  can  be  determined  if  the  density  n(r^) 
at  a  radius  r^  is  known.  Then 

A     =  r       ,   — -  ,"(r' -   w  (3-72c) 

The  radial  current  density  is 

rr^LAj£L  {  J,  (y£r)4-  CM  (JZJ]  (3-3) 

d.  Coaxial  Cylinders 

The  equations  of  the  previous  section  also  describe  the 

particle  and  current  density  relations  which  exist  between  the  surfaces 

of  coaxial  cylinders  of  infinite  length  with  inner  radius  R  and  outer 

radius  R1  .  Because  of  the  boundary  conditions,  )/   is  no  longer  an 

independent  variable,  but  is  an  eigenvalue  in  the  steady  state  plasma. 
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Using  the  boundary  condition  of  zero  density  at  the  walls,  the  value 
of  JJL    is  seen  from  equation  (3=72)  to  be  that  value  for  which 

£■(£/)  +  C  Nol&Jt)  =  O  (3-7W 

h»     Spherical  Geometries 

Distribution  functions  for  decaying  plasmas  with  arbitrary- 
initial  distributions  and  for  steady-state  plasmas   in  spherical  geo- 
metries can  be  obtained  by  series  development  of  equations   (3-17). 
Equations  for  the  spherically  symmetric  plasma  in  the  steady  state 
(with  ionization  occuring)  are  very  similar  to  (3-66)  and  (3-72)   with 
spherical  Bessel  functions  (equation  (3-17c))   replacing  the  cylindri- 
cal Bessel  functions   in  equations  describing  spherical  geometries  corres* 
ponding  to  the  cylindrical  geometries  treated.     In  addition,    it 
is  of  interest  to  describe  the  steady  state  solution  for  a  steady  or 
very  slowly  decaying  plasma  in  which  there  is  no  ionization  occuring, 
and  in  which  a  sphere  has  been  placed.     We  treat  the  latter  problem 
first. 

a.     Sphere  in  Homogeneous  Plasma-no   ionization  occuring 
If  we  have  a  spherically  symmetric   infinite  plasma 
with  an  undisturbed  particle  density,   n_,    in  which     l) =  $&    —  q  , 
with  neither  magnetic  nor  electric  fields    imposed,  we  have  from 
equations   (3-17) (with  k/>  -  m  -  0)   the  following  density  distribution 
around  a  sphere  of  radius     R 


7t(r)*  y)0-  7T  r^  (3~75) 
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We  have  employed  the  boundary  conditions  that  n(r)  be  zero  at  the 
sphere  and  equal  n0  at  large  distances.     The  zones  of  influence 
for  such  a  sphere  are 

90%    -      r±\oRt  (3-76) 

and  the  radial  current  density  is 

P   -   —  Q> ytc  ^  (3-77) 


r  - 


r 


2. 


/  noDq 

This  current  density  is  a  maximum  (of  magnitude     )  at  the  surface 

R 

of  the  sphere  and  approaches  zero  as  l/r^  at  large  r.  The  current  is 
directed  toward  the  sphere  and  represents  a  flow  of  particles  from  a 
source  at  infinity  to  a  sink  at  the  sphere. 

The  heat  liberated  at  the  surface  of  the  sphere  is 

£«Ut"    fisR^Dj^  +  f-).  (3-78) 

It  is  proportional  to  the  square  of  the  radius  of  the  sphere. 

Using  the  equations  thus  stated,  it  should  be  possible  to  calculate 
the  particle  density  distribution  which  a  small  spherical  surface  in 
a  large,  slowly-decaying  plasma  would  establish  about  itself, 
b.  Spherical  Shell~internal  plasma 

From  equations  (3=17)  it  is  also  possible  to  obtain 
easily  the  steady  state  relations  for  a  plasma  bounded  by  a  sphere  of 
radius  R  centered  within  a  spherical  shell  of  radius  R1.  We  shall 
again  assume  spherical  symmetry *  thus  limiting  ourselves  to  the 
spherical  Bessel  functions  of  zero  order.  By  employing  the  follow 
ing  relations  for  these  functions 

c^(kr)  (3-79) 


Yto(kr)    -    - 
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and  doing  some  manipulating*  we  obtain 

W(r)   =  C,  r'W^LIVr^'J     %'>r>K)  (3-80) 

vhere 

%  --  h  s(t&£  (3-80a) 

and  for  which  derivation  the  condition  of  vanishing  density  at  the  sur- 
faces was  employed.  Ci  is  an  arbitrary  normalization  constant  which 
can  be  determined  if  the  density  is  known  at  some  radius  r^  within 
the  plasma. 
Then, 

and 

C-     %n    ,sc[^|].  (3-80O 

The  radius  of  maximum  density  R» * ,  can  be  found  by  setting  the  first 
derivative  of  (3~80)  to  zero,  vhich  yields  the  following  equation  for 
R»« 


•*WSf*l  •  Wttm 


R^i? — Lfe'-Bj  J  "   *"  "  !-«-•«  j.         (3.81) 


For  a  sphere  within  an  infinite  plasma  in  which  ionization  is  occuring* 
this  is  the  radius  at  which  the  effect  of  the  spherical  surface  be- 
comes zero. 

The  current  density  within  this  plasma  is 

*  \FUFr)c<t*Ir' -  ><  -1        r         L^r^J      (3.82) 

It  is  to  be  noted  that  it  has  two  components.  The  first  is  directed 

toward  both  surfaces  from  a  radius  within  the  plasma,  and  the  second 
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is  directed  radially  outward  from  the  center,  but  is  zero  at  both 
surfaces. 

For  a  spherical  shell  of  radius  R'  vith  no  internal  sphere,  one 
obtains 


7f(r)   -  ^T   3^(-g7J 


(3-83) 


with 


J-       JL      JL 

A^  Ds  fc'  (3-83a) 

This  equation  was  obtained  from  equation  (3-80)  by  setting  R  ■  0}    and 
by  evaluation  of  C,   by  L'Hospital*s  rule 


(3-83b) 


r+o   p r-*o <-,  R/  s  r\o 

where  nQ  is  the  density  at  the  center.  Figure  3-6  is  a  graph  of  this 
function. 

5.  Prolate  Spheroidal  Geometries  (the  wire  or  probe) 

The  influence  of  a  prolate  ellipsoid  upon  a  plasma  is  nov 
considered  (see  Appendix  I  for  transformation  equations  and  figure). 
In  its  limit5  when  the  minor  axis  is  small  and  the  major  axis  approaches 
the  length  d,  the  effects  of  such  a  body  can  be  used  as  an  approxi- 
mation of  the  effects  of  a  finite  wire  or  narrow  probe  of  length  d. 

Flammer  [6]  shows  that  the  wave  equation 

(72+  k2)  N  =  0,  ( 3-81*) 

equation  (3-7)5  expands  in  prolate  spheroidal  coordinates  toj 

62 


O  .2.  W  .(<>  <%  \.0 

Radial     ?OS\hon     (&<xch'cm  cS  R") 


flhure    3~6>       Rocliai   D<?nsif</   l>sfnW*;!rrA  {or  Vta|bLO$pk*ra 


63 


in  which  we  have  set 

C-ikd;  (3-8$a) 

and  that5  by  the  method  of  separation  of  variables,  solutions  to 
(3-8£)  may  be  obtained  in  the  form  of  the  Lame  products 

Ak„  «    S*&f)  £«.  tCj  t^W  (3-86) 

The  functions  Smn(c/?7)  and  Rmn  (c,£)   satisfy  the  ordinary  differential 

JGj  [O-D&i  ^(c>H  4A-"c'f  -  7f\SmAj>  -  O     (3.86a) 

^k?H)^B(^)]-[A--'Y  +^f]^?)=°  (>86b) 

in  which  the  separation  constants  <\mn  and  ft\   are  the  same  in  both 
equationso  It  is  a  combination  of  such  functions  (of  the  first  and 
second  kinds)  that  is  given  in  Equation  (3-18)  as  the  general  solu- 
tion to  (3-81*) 

a.  Prolate  Ellipsoid  in  stationary  infinite  plasma 

Before  proceeding  further  with  the  prolate  wave  func- 
tions, one  special  case5  analogous  to  that  of  the  sphere  in  the 
infinite  stationary  plasma,  should  be  treated.  By  imposing  the 
boundary  condition  that  at  large  £  the  density  n  of  a  stationary, 
non-ionizing  plasma  be  independent  of  l^  and  7f  (so  that  §**-  -  0 
and  -Tzpi.  ~  0),  we  can  reduce  equation  (3-85)  to 


^[cf-o^]  -o 


which  integrates  to; 


(3-87) 


VI-  >l\\    -    C.°^'^\\  i<*?<LV<00  (3-88) 
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vhere  n  is  the  uniform  density  at  infinity,  and  the  density  at  the 
o 

surface  T  is  zero0  The  current  density  is  then? 

n  ^ 

To  investigate  the  properties  of  this  solution,  vre  shall  use 
the  identity 

eMTXCi  =  i^-n)  -  ilU(]-l)  (3-90) 

by  which  Equation  (3=88)  becomes? 


71 '-*L     it~cu\)-±*(i>-i) 


(3-91) 


at  very  large  distances  (i.e.,   <T  >>  1),  In  (T  +  1)  -  ln(?  -1)  ~  -r 
and  equation  (3-91)  can  be  expanded  to 

71  *   *e0  *  ^)  (3-92) 

which  is  identical  to  equation  (3-75)  for  spherical  geometry,  a 
result  which  reflects  the  fact  that  the  surface  "looks  like"  a  sphere 
at  large  distances.  In  the  vicinity  of  the  surface,  i.e.,  for  ^^^e, 
we  let 

(with  £  such  that  the  mean  free  path,  A  ,  is  much  less  than  €   ), 
and  expand,  obtaining 
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*ft 


(3-9lib) 


&-9Uc) 


By  neglecting  second  order  terms  in  £  and  using  the  expansion 
equation  (3-9lic)  can  be  reduced  to 


(3-95) 


7i  #  -yio 


?.H  V 


-  Tfc 


£€ 


(3-96) 


Thus,  for  the  density  at  the  surface,  ^  -  0  and  therefore  n  =  0.  Equa- 
tion (3-96)  shows  that  the  density  is  a  linear  function  of  the 
distance  fron  the  surface  for  short  distances. 

At  large  distances,  the  current  density  can  be  seen  from  Equa- 
tion (3-92)  to  be  an  inverse  square  function  of  distance  (as  it  is 
for  a  sphere)? 

(3-97) 


At  short  distances,  we  obtain  from  (3-96) 


r 


(3-98) 

which  is  a  minimum  at  the  center  of  the  surface  (y)  ■  0)  and  approaches 

infinity  at  the  ends  (¥1   ■  £  1).  This  occurs  because  the  area  at  the 

ends  is  zero,  and  thus  any  current  flowing  towards  them  constitutes 

£6 


an  f infinite*  current  density. 

The  current  density  at  the  surface  is  found  front  (3-98)  to  be 

which  again  becomes  infinite  at  the  ends  and  is  a  minimum  at  the 
center.  Ke  can  obtain  an  estimate  of  the  heating  at  the  surface 
per  unit  length  if  we  approximate  the  ellipsoid  with  a  cylinder  of 
radius  r  ■  (  \    -  1)  and  approximate  the  current  density  with  that  at 
the  center  (/I  -   0). 

£^1    >_  *   <*    *,*.(•*+%-)  (3.100) 

To  find  the  behavior  at  small  radii,  we   let     5  «*  1  +  £ft     obtaining; 

which  approaches  zero  as  £  approaches  zero. 

b.  Prolate  Spheroid  in  steady  state  plasma  -  ionization 
occur ing 

Given  an  infinite  plasma  of  homogeneous  density  n0  at 
large  distances  and  a  prolate  spheroid  whose  surface  is  the  surface 
£,   ,  one  can  simplify  equation  (3— 18)  considerably  by  consideration 
of  boundary  conditions.  First,  the  requirement  that  the  wave  function 
be  finite  at  71   «  +  1  confines  the  solution  to  angle  functions  of 
the  first  kind  since  the  angle  functions  of  the  second  kind  are  singular 
at  these  points.  If  we  look  for  those  solutions  for  which  the  den- 
sity, n;is  uniform  at  large  distances  for  all  f  ,  then  N  is  a  function  of 

>l   and  ^  alone,  and  m  ■  0.  Further,  the  prolate  angle  functions 
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reduce  to  the  associated  Legendre  functions  when  c  vanishes;  i.e., 
as  the  spheroidal  insulator  becomes  spherical  in  shape,  Therefore, 
we  can  determine  the  subscript  "n"  in  Equation  ( 3 - 1 3 )  by  studying  the 
Legendre  function,  P^(  f)  )   because  our  solution  must  satisfy  the  limit- 
inn  case  of  a  hemispherical  body.  Only  Pq(  77  )  is  positive  for  all  7\ 
in  the  range  -1^  7\  <   1.  Because  n  must  be  positive  for  all  >)  ,  the 
subscript  must  be  0.  Consequently,  we  obtain  for  the  density  distri- 
bution 


n  =  n.  Si?  (crflBoa  (c, f )  1 3 R»  (.c,  (i  j,  ( 3-102) 


where 

r        ^i         dlZ  (3-102a) 

If  we  employ  boundary  condition  (2-99)  (^vall  -   °)> 

b  =  -  RSctf)  °~102b) 

It  must  be  remarked  that  this  solution  does  not  satisfy  the  require- 
ment that  the  density  be  independent  of  ))   at  large  distances*  However, 
for  small  values  of  c,  and  therefore  for  either  very  small  or  for  nearly 
spherical  bodies  the  71  dependence  is  negligible.  The  deviation  is 
only  )<%   for  c  ■  0,5,  but  increases  to  practically  95%   for  c  -  5,0. 
c »  Prolate  spheroid  and  infinite  va 1 1 

Let  us  consider  a  semi-infinite  plasma  of  uniform  density 
n  at  large  distances,  bounded  only  by  a  single  infinite  plane  which,  in 
prolate  spheroidal  coordinates,  is  Y)  =  0.  Further,  consider  an 
insulator  (probe)  whose  surface  is  defined  by  "<-  Kb  pro jec ting  into  the 

plasma  as  indicated  in  Figure  (3-7).  The  solution  to  the  continuity 

68 


tl 


Fmure    3'?      Prolate    Spheroid    an  A  InffWik  UloJ 


69 


equation  is  developed  in  the  sane  manner  as  in  the  previous  section, 
except  that  in  this  instance  the  appropriate  Legendre  function  is  Pj|(  Y(  ) 
because  it  has  the  property  of  being  positive  for  all  values  of 
in  the  range  0  to  +1  and  further  of  being  zero  at  71 «  0.  Thus  the 
density  is 

7}^7?oSc,(c^)[Ro,(c,\)  4    BR^fe/^J^  (3^03) 

where,  again, 

/»=  y  s    4/?  (3-103a) 

and  if  %all  s  °> 

K0.fc)j»)  (3-103b) 

b  "  "   £o>^ 
The  field  and  potential  equations  along  the  axis  of  the  probe 

(i.e.,  for  yi*>  1)  ares 


/p  \         ,  z  D+-h_     LKo,  ^^-h^rso,  cc^jj 

ff<  ^  " J  ^  [r^)-b^)J  °   4) 

where  the  prime  indicates  the  first  derivative,  and 

^v'2.;?j*/   T  /a+-w-^~  *,  (3-io5) 

6.  Oblate  Spheroidal  Geometries  (circular  lamina) 

In  general,  the  solutions  to  the  continuity  equation  in 

oblate  spheroidal  coordinates  are  the  same  as  for  the  prolate  with 

-ic  replacing  c  and  +i£  replacing  5  in  the  equations.  Consequently 

these  will  not  be  further  discussed,  except  for  the  problem  of  the 

circular  lamina  in  an  infinity  stationary,  non-ionizing  plasma. 
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a.  Lamina  in  stationary  plasma 

To  study  the  influence  in  an  infinite  plasma  of  a 
circular  lamina,  the  effect  of  an  oblate  ellipsoid  is  first  discussed 
and  then  the  limit  as  the  semi-minor  axis  approaches  zero. 

Proceeding  as  for  the  prolate  case,  by  imposing  the  condition  that 
at  large  £  the  plasma  density  be  independent  of  if  and  yj  ,  and  by  sub- 
stituting -H  ?  for  V  ve  obtain  from  equation  (3-85) 

%[<»?)%]-<>•  (3-06) 


Integrating  ti-dce $  ve  obtain 


?l  =  c,  w'Y?)  +c^ 


(3-107) 


Letting       -£^  97   =  fc   ,  we  get 

7,0  --  C,  h^'(co)  -fCa  =  hSV  C,-*C,_  (3-ioB) 

Using  as  the  boundary  condition 

%(?>)  =  0  =  C,  -hvZ'Cft  +Ct  (3.109) 

^subtracting  this  from  equation  (3=108)  yields 

so  that,  for  an  oblate  ellipsoid  in  an  infinite  plasma, 
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VJhen  the  semi-minor  axis  is  zero,     T  -  0,  and  the  figure  is  a 
circular  lamina  of  radius  rQ  -  c.     Equation  (3-12£)  becomes 

"*         I.S7\~~~  (3-U3) 

The  zone  of  influence  of  this     amina,  within  vhich  n  ^••J-nQ^is  given 

by 

i.e.,  "jT^  (3-llteO 

which  is  an  oblate  ellipsoid  with  semi-major  axes  r  ■    N/2r0  and 
semi-minor  axis  z  =  rQ0 

The  current  density  in  this  plasma  (  with  c  ^  fQ   )  will  bei 

This  current  density  is  a  maximum  at  the  surface  of  the   lamina, 
decreases  as     -r,  to  zero  towards  infinity,  and  because  the  area  of 
the  "side"  of  tte    lamina  is  zero,   is  infinite  at  the  sides  (  71  =0)  of 
the  "surface".     If  we  assume  the  current  density  at  the  center  of  the 
surface  to  be  an  approximation  to  the  surface  current  density  across 
the  entire  surface,  then  the  recombination  energy  liberated  per  unit 
surface  will  be  approximately 


c  ~    7Tlkyiafo  /  w    .    leT_\ 

£5^€  *  Tst7 —  (eV{  f  "S"  /  ( 3-ll6) 

which,  as  for  the  sphere,  is  proportional  to  the  radius. 
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IV 
APPROXIMATE  SOLUTIONS  FOR  INSULATED  PROBES 

Mien  an  insulated  probe  is  inserted  into  a  dense  plasma— that 
is,  a  plasma  in  vhich  the  mean  free  path  is  equal  to  or  less  than 
the  probe  dimensions— Langmuir' s  probe  theory  £8]  may  not  be  appli- 
cable. Particle  losses  at  the  probe  surfaces  may  cause  appreciable 
effects  on  the  particle  density  distribution,  and  hence  on  the  probe 
potential.  Under  these  conditions,  o£   of  equation  (2-71)  can  be 
appreciably  large  and  consequently  so  may  the  effect  of  the  in- 
sulator on  the  steady-state  distribution.  The  approximation  describ- 
ed in  Section  III-A~2  has  been  employed  to  calculate  this  effect 
for  a  floating  probe  in  a  partially  ionized  gas.  Tvo  different 
geometries  are  considered.  The  first  is  a  prolate  llipsoid 
inserted  radially  into  a  cylindrical  discharge,  vith  the  sensing 
element  on  its  tip.  The  other  is  a  cylindrical  probe  aligned  parallel 
to  the  discharge  axis.  For  the  latter,  ve  are  able  to  account  for 
the  effects  of  a  longitudinal  magnetic  field. 
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A.  The  Prolate  Probe 

We  do  not  have  an  exact  solution  of  the  continuity  equation 
(2-55)  for  a  cylindrical  or  ellipsoidal  probe  projecting  from  a 
cylinder  wall  into  a  plasma.  However,  we  have  seen  that  there 
is  a  solution  for  the  density  distribution  in  a  long,  narrow 
plasma  (equation  (3-66),  In  the  absence  of  a  magnetic  field,  it  is 


n-n0^o(^)) 


(U-l) 


in  which 


(U-la) 


—  H 


E 


3Ure  *H        QncAate    \>ra\>e     m     C^lwdrica]    Discharge 
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Further,  we  have  obtained  a  solution  for  a  prolate  spheroid  attached 
to  an  infinite  wall  (equation  (3-102)) 

*  *  «,'  So,  (C, f )  [ffife,  f )  i-  B  Ko?&,  I)]     ( u-2) 

Based  on  the  assumptions  outlined  in  Section  III-A-2,  we  shall 
approximate  the  effect  of  a  narrow  cylindrical  or  prolate  ellipsoidal 
probe  by  joining  the  two  models  described  above  as  indicated  in 
Figure  h-1.  The  infinite  wall  is  tangent  to  the  cylinder,  and  the 
z-axis  of  the  probe  (  7j   ■  1)  lies  along  the  diameter  of  the  cylinder. 
We  assume  that  the  presence  of  the  probe  does  not  significantly  alter 
the  ionization  or  diffusion  rates,  so  that  equation  (3-103a)  becomes 

where   —  jsii  ,  the  length  of  the  probe  (if  the  minor  axis  is  much 
shorter  than  the  major  axis)  and  R  is  the  radius  of  the  cylinder. 
No  effort  will  be  made  to  describe  the  density  throughout  the 
cylinder.  We  shall  merely  attempt  to  predict  the  density  distri- 
bution along  the  7)  s  1  axis  of  the  probe  (i.e.,  along  the  diameter 
of  the  cylinder).  It  is  this  density  distribution  which  is  res- 
ponsible for  the  potential  distribution  sampled  by  a  sensing  element 
at  the  tip  of  the  probe.  In  accordance  with  the  procedure  outlined 
in  Section  III-A-2,  we  readily  obtain  the  following  condition  for 
determining  the  position  along  this  axis  of  the  dividing  surface  be- 
tween the  two  models 

\ta.\  -1^1 

1  ^  'Probe     '"*  Cylinder  (h-3) 

Equation  (h-3)  was  solved  graphically  because  sufficient  published 

data  on  the  derivatives  of  the  spheroidal  wave  functions  was  not 

available. 
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Because  of  this  paucity  of  information  on  the  spheroidal  vave 
functions,  it  was  also  necessary  to  limit  the  study  of  this  problem 
to  the  limiting  case  where  A  is  so  small  that  the  boundary  con- 
dition (2=99)  (n  ..  -  0)  may  be  employed.  For  this  problem,  B  of 
vail 

equation  (h-2)  is  readily  obtained! 

The  normalization  constant,  n',  of  equation  ( U— 2)  was  determined 

o 

in   relation  to  n  of  equation  (h-1)  by  the  graphical  solution  of 

equation  (U-3) • 

A  study  of  the  solutions,  and  of  the  physical  properties  of  our 

model,  reveals  that  this  model  can  be  justified  as  an  approximation 

to  the  particle  distribution  in  the  vicinity  of  the  probe  for  probe 

penetrations  varying  from  zero  to  about  half  the  cylinder  diameter. 

For  short  probe  lengths,  particle  losses  to  the  adjacent  regions  of 

the  infinite  vail  (which,  for  the  short  probe,  is  a  reasonable  approxi- 

mation  of  the  adjacent  cylinder  wall)  is  the  predominant  factor  in 

establishing  the  particle  distribution  in  the  vicinity  of  the  probe. 

For  instance,  for  a  probe  length  of  0.2R,  the  parameter  c  of  equations 

(U-2)  is  0.1$.  From  Figure  11=1  (Appendix  II)  we  see  that  the  pro- 

late  angle  function  S^''  is  very  nearly  that  for  c  =  0,  which 

corresponds  to  zero  probe  penetration.  Consequently,  the  effect  of 

the  short  probe  is  very  slight  in  comparison  to  the  effect  of  the 

wall,  as  we  expect  from  consideration  of  the  surface  which  would  be 

the  cause  of  most  of  the  particle  losses  from  the  region  around  the 

probe.  However s  for  deeper  penetration,  the  angle  function  changes 

appreciably,  indicating  that  the  surface  of  the  probe  exerts  a  more 
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dominant  influence.  Again,  this  is  expected  since  particles  in  the 
vicinity  of  the  tip  of  the  probe  are  more  likely  to  diffuse  to  some 
portion  of  the  probe  than  to  succeed  in  reaching  the  vail.  Thus, 
even  though  the  infinite  vail  is  no  longer  a  good  approximation  of 
the  cylinder  wall  for  the  deeper  penetration,  the  increased  probe 
surface  losses  justify  using  the  resulting  particle  density  distribu- 
tion functions  to  describe  the  plasma  near  the  end  of  the  probe. 
When  the  probe  penetrates  more  than  half-vay  through  the  plasma, 
the  influence  of  the  opposite  vails  of  the  cylinder  on  the  particle 
distribution  near  the  end  of  the  probe  become  increasingly  more  important, 
and  our  model  breaks  dowi. 

Conversion  from  prolate  coordinates  to  probe  dimensions  is  accomplished 
vith  the  fol loving  equations. 

probe  length?  (  =  ~g  f0  ~'  (U-S) 

probe  radius  (  =  semi-minor  axis):  rp  =  ■•> J^- 1  -  ^r J r L T *f /£- 1( I1-6 ) 
^Tr  X      is  the  equation  of  the  probe  surface. 

Flammer  [6~]  vas  the  source  of  the  information  on  the  spheroidal 
vave  functionSo  Data  not  calculated  or  graphed  in  this  reference  is 
presented  in  Appendix  II. 

The  graphs  on  the  following  pages  vere  computed  using  the  proce- 
dures outlined  in  this  section  and  subsidiary  equations  developed  in 
preceeding  sections. 

Figure  U-2  is  a  graph  of  density  versus  distance  for  probes  of 
several  radii  but  equal  depths  of  penetration.  Note  that  the  magnitude 
of  the  dependence  of  the  probe  influence  on  the  probe  radius  is 
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relatively  small,  an  order  of  magnitude  change  in  probe  radius 
causing  an  increase  in  the  region  of  influence  of  only  £C$. 

Figure  Ii~3  is  a  graph  of  the  density  distribution  for  various 
depths  of  penetration  of  a  probe  of  radius  r_  =  0.02R  (which,  for  a 
2.5  cm  discharge,  corresponds  to  a  0.5  mm  probe). 

Figure  h-k  shows  "effective  probe  length"  vs.  true  probe  length 
for  the  0.02R  probe.  The  effective  length  is  based  on  the  inter- 
section of  the  two  models,  and  is  approximately  the  point  of  the  ori- 
ginal particle  distribution  that  the  probe  samples.  The  curve  marked 
X  *>f      2r  is  the  effective  probe  length  when  the  probe  does  not 
disturb  the  plasma,  and  corresponds  to  the  true  probe  length.  The 
broken  segment  of  the  curve  is  an  estimate  for  that  region  in  which 
this  model  breaks  down,  as  outlined  above. 

Figure  k-5   is  a  graph  of  the  potential  disturbance  associated 
with  various  depths  of  probe  penetration.  This  potential  disturbance 
is  defined  as  the  potential  drop  between  the  point  the  plasma  samples 
(the  effective  probe  length)  and  the  sheath  edge.  A  sheath  depth  of 
0.01R  was  assumed  for  this  calculation.  These  potentials  can  be 
used  with  the  effective  length  data  to  analyze  experimental  results 
or  to  predict  them  from  the  theory. 

Figure  l*-6  represents  such  a  potential  prediction.  It  is  the 
prediction  of  the  potential  measured  by  such  a  probe  assuming  that  it 
actually  measures  a  difference  between  the  potential  at  the  maximum 
density  i eg ion  of  the  plasma  (zero  potential)  and  the  potential  at 
the  sheath  edge  at  the  tip  of  the  probe. 

A  very  significant  prediction  of  the  above  theory  is  that  the  tse 
of  an  insulated  probe  in  a  plasma  of  short  mean  free  path  will  cause 
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a  definite  asymmetry  even  in  an  originally  symmetric  plasma.  Although 
these  results  are  not  applicable  directly  to  a  plasma  in  the  presence 
of  a  magnetic  field,  some  of  this  effect  is  still  expected. 
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B.  Cylindrical  Probe  Parallel  bo  Axis 

Figure  li-7  is  a  top  view  of  a  possible  insulated  probe  con- 
figuration for  measurements  in  the  presence  of  a  magnetic  field.  The 
cylindrical  probe  (enlarged  in  the  figure)  is  supported  at  its  ends 
and  long  enough  that  the  effect  of  the  end  supports  are  negligible 
in  the  region  at  the  center  of  the  probe.  The  sensing  element  (or 
elements)  are  located  at  the  center  of  the  probe  on  the  common 
diameter  of  the  probe  and  plasma. 

There  is  only  one  alignment  of  the  probe  for  which  an  exact 
solutionaf  the  continuity  equation  (2-5>7)  is  possible:  when  probe 
and  cylinder  are  co-axial.  As  previously  developed  (equation 
(3-72)),  the  solution  is 

■n-*'{tilij\  +  <^<£/>]        rt±r<R       (u-7) 

with  r  the  probe  radius,  R  the  cylinder  radius,  and  n1  the 
normalization  constant.  If  A  is  not  very  much  less  than  the  probe 
dimensions,  boundary  condition  (2-71)  should  be  used  at  the  probe 
surface.  However,  equation  (2-99)  (n  i,  *  0)  is  often  still  appli- 

Wet  I  Jl 

cable  at  the  cylinder  wall,  and  was  used  for  our  calculations. 
Using  equation  (2=71) $  which  is 


at  the  sheath  yields,   from  (li-l) 


(U-8) 


(k-9) 


8U 


Figure     4-7      CyUplrltAl    fi-<\>e    PavcxWA    \n    <',/ 
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from  which  we  obtain 

Next,  employing  equation  {2-99)   at  the  outer  cylinder,  we  obtain 

£(§£*)  =  CN0(J&J<)  ih'n) 

Equations  (U=10)  and  (U=ll)  are  sufficient  to  determine  both  C  and 
Jy/r-,  and  they  are  most  easily  solved  by  successive  approximation. 

Table  L-l  lists  the  i  nese  parameters  for  various  p£ between 

zero  and  one,  using  the  following  probe  dimensions  and  mean  free  path 

r  •  0,02R 
P 

X    s  0.02R  (h-12) 

(again  this  corresponds  tp  a  0„5  mm  probe  in  a  2.5>  cm  cylinder). 
Even  when  there  is  no  recombination  at  the  probe  (  C^.   ~   0),  the  rate 
of  ionization  for  the, steady  state  must  be  greater  than  when  the  probe 
is  not  present  because  the  probe  still  occupies  a  finite  vo  lume  in 
the  cylinder,, 

When  the  .probe  is  off~axiss  it  is  necessary  to  employ  the 
approximation  described  in  Section  III-A-2.  Equation  (h-3)  is  still 
applicable  for  determining  the  separation  of  the  zones  along  the 
common  diameter  (the  only  location  in  which  we  will  be  interested). 
In  order  to  evaluate   J'%^  *n  the  vicinity  of  Jthe  probe  for  various 
probe  positions.,  we  interpolate  between  the  calculated  value  when  the 
probe  is  coaxial  and  the  value  it  would  have  when  the  probe  is  at 
the  cylinder  wall  (and  has  negligible  effect  on  the  eigenvalue  because 
of  the  low  density  there).  This  latter  value  we  assume  to  be  th«~ 
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familiar  2„ii05/R»  Assuming  that  this  parameter  would  depend  on  the 
density  of  the  plasma  in  the  position  in  vhich  it  is  placed,  we 
obtain  the  following  relationship 


T         I  r         ~n(r) 


(U-13) 


where 


From  the  values  thus  calculated^  a  mean  was  calculated  for  each  &C 
(see  Table  U-l)*  and  this  value  was  used  in   the  computations  of  the 
density  distribution  in  the  vicinity  of  the  probe.  The  use  of  a 
more  precise  value  is  not  justified  by  the  nature  of  the  approximation. 
For  the  density  distribution  in  the  regions  controlled  by  the  cylinder 
walls,  the  value  of  J^/r,        for  the  cylinder  alone  was  used  (2,ho5/R) 
because  the  exact  effect  of  the  probe  on  this  distribution  could  not 
be  accurately  estimated.  The  degree  of  error  resulting  from  these 
approximations  can  be  inferred  from  a  comparison  of  the  exact  solution 
for  the  coaxial  probe  and  the  approximate  solution  shown  in  Figure 
U-llj.  As  expected  from  the  approximation  model,  the  maximum  error  is 
in  the  region  where  the  two  models  are  joined  (where  the  positional 
error  is  0.1R  and  the  magnitude  of  the  approximation  is  12%   greater 
than  that  of  the  exact  solution) .   The  magnitude  of  the  error  de- 
creases towards  either  surface,  but  the  percent  error  remains  approxi- 
mately constant.  Consequently,  the  potential  determinations  (which 
are  based  on  a  logarithm  of  a  ratio  of  densities)  will  be  reasonable 
approx  ima  t  ions « 
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Several  distribution  carves  have  been  obtained  by  using  the 
procedure  outlined  aboves     (1)  Figure  li-8  is     a  graph  of  the  exact 
solution  for  a  0.02R9   co~axial  cylinder  for  various   U.  j   (2)Figure 
h-9   is  a  graph  of  the  exact  solution  for  co-axial  cylinders  of  various 
radii  at  a  given  o<j   (3)  Figure  li-10  is  the  density  distribution  for 
various  ££  for  a  0.02R  probe  in  an  off -axis  posit ion^   and  (U)  Figure 
h-11  is  for  a  given  c<  and  various  probe  positions  for  a  0.02R8  probe. 

Effective  Probe  Position  vs.  actual  probe  position  is   indicated 
in  Figure  h-12  and  potential  correction  to  be  applied  to  probe  readings 
in  Figure  li-13j  we  have  assumed  the  sensing  point  to  be  on  the  side 
facing  the  center  cf  the  plasma  when  the  actual  position  is  less 
than  }R  from  the  wall.     Potential  curves  are  given  in  Figures  It-lli 
and  h-l5«     The  former  one  is  for  one  contact,,  oriented  as  indicated 
above.     The  latter  is  for  two  contacts,  one  on  either  side  of  the 
probe  along  the  common  diameter.     Figure  ii-13,  the  potential  correction 
curve,  is  just  the  potential  difference  between  the  sheath  edge  and  the 
point  of  intersection  of  the  zones  of   influence  for  various  cC  and 
probe  positions. 

Figure  1|-16  is  a  curve  of  heat   liberated  at  the  probe  for 
various  values  of  probe  radius  when  the  probe  is  aligned  with  the 
center  of  the  cylinder,  and  li=17  is  for  various  values  of      oL   .     It  is 
assumed  for  these  calculations  that  the  fraction  of  particles  which 
diffuse  to  the  sheath  reach  the  surface  and  recombine,   liberating 

(eVc  ■*■  kT%  )  of  enersy° 

The  results  of  the  above  theory  indicate  that  the  cylindrical 
probe  should  have  considerably  less  effect  on  the  particle  distribu- 
tion functions  than  the  prolate  probe  in  a  given  cross-section  of 

the  plasma.  This  would  be  expected  from  the  decreased  surface 
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involved  for  the  cylindrical  geometry „  Consequently,  there  is  less 
assymetry  produced  by  the  probe;,  and  the  potential  curves,  particularly 
vith  contacts  on  both  sides  of  the  cylindrical  probe,  should  be  more 
symmetric*  However 5  the  effect  of  measuring  a  potential  at  an 
"effective  distance'8  causes  some  flattening  of  the  predicted  curve. 
Figjre  h-9   indicates  that  the  density  distribution  in  the  vicinity 
of  the  probe  is  relatively  insensitive  to  the  probe  dimensions. 
Consequently,  the  potential  drop  associated  vith  this  probe  configura- 
tion is  relatively  constant  for  various  probe  radii. 
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TABLE  U-l     CYLIMDRKAL  PROBE  PARAMETERS— EXACT  VALUES  OF 


J~^n?  and  C       Vs.     oi.      for  probe  at  center,  and  mean 

values  vs.    oC, 


oL 

Probe  at 

Center 

Mean  Values 

np 

C 

n/d^j. 

C 

i 

2.91 
R 

.569 

2.70 
R 

.5514 

»5 

2o8h 
R 

ol4?8 

2,66 

R 

.U67 

.25 

2.75 

R 

.365 

2.61 

R 

.361 

.1 

2,6l 
R 

,2l6 

2.53 
R 

.215 

.05 

2.53 

R 

.131 

2.1|8 

R 

.131 

.02 

2j£? 
R 

.06h 

2.U4 
R 

.063 

,01 

2.J4h 

R 

,037 

2.U3 

R 

.036 

0 

2oiil 

R 

.00? 

2.U1 

R 

.007 
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Co  Evaluation  of  Prohe  Measurements 

In  a  dense  plasma,  in  t-hich  the  mean  free  path  is  of  the 
order  of  magnitude  of  or  less  than  the  probe  dimensions^  strong 
deviations  from  Langmuir»s  theory  can  be  expected.  These  deviations 
depend  on  the  size  and  character  of  the  probe*  By  using  the 
approximation  described  in  this  thesis,  one  can  calculate  these  effects 
in  terms  of  an  "effective  probe  length'1  and  a  "potential  disturbance^ 
These  two  parameters  can  then  be  used  to  analyze  or  predict  the 
experimental  results  in  the  manner  just  employed  for  the  prolate  and 
cylindrical  probes. 
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V 
MAGNETIC  FIELD  EM)  EFFECTS 

Let  us  now  investigate  the  effects  on  the  particle  density 
distribution  functions  of  the  decrease  in  magnetic  field  intensity 
st  the  end  of  a  coil  with  attencent  intersection  of  field  lines  with 
the  container  walls.  The  effect  expected  is  that  diffusion  along  the 
field  lines,  being  greater  than  that  perpendicular  to  the  vails  (radial 
diffusion),  would  cause  the  density  distribution  near  the  ends  of  the 
field  coi)  to  decrease  more  with  increasing  radial  distances  from  the 
axis  than  predicted  from  radial  diffusion  alone.  Particle  losses  by 
longitudinal  diffusion  to  the  end  of  the  field  coil,  where  some  of 
the  field  lines  intersect  the  container  vails,  would  be  responsible 
for  this  effect,  which  would  manifest  itself  physically  as  a  con- 
striction near  the  ends  of  the  field  coil. 

An  exact  solution  of  the  continuity  equation  for  this  configuration 
(Figure  5-1— note  the  enlarged  radial  scale)  could  not  be  found,  and 
consequently  the  approximation  described  in  Section  III~A<=2  was  used. 
The  plasma  is  divided  into  zones  in  which  the  diffusion  would  be 
primarily  radial  or  primarily  longitudinal  (along  field  lines). 
The  latter  zone,  in  the  vicinity  of  the  intersection  of  the  field 
lines  with  the  plasma  container,  is  a  zone  in  which  the  particle  dis- 
tribution is  controlled  by  the  loss  of  particles  diffusing  along  the 
field  lines  and  recombining  at  the  region  of  the  wall  where  the  lines 
intersect  the  wall.  The  size  of  the  region  of  parallel  diffusion 
increases,  of  course,  with  increasing  magnetic  field.  If  we  neglect 
the  curvature  of  the  fielc  lines  and  use  the  boundary  condition 
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we  obtain  for  the  density  distribution  along  a  field  line 

where  n"  is  a  normalizing  factor  to  be  detetminci,  the  z-axis  is 
the  cylinder  and  magnetic  field  axis,  and  the  origin  is  the  inter™ 
section  of  the  field  line  and  the  wall.  This  solution  cannot  be 
rigorously  derived,  but  is  justified  as  an  approximation  consistent 
with  the  accuracy  of  this  calculation.  Comparison  with  the  solution 
for  a  bounded  cylinder  in   a  homogeneous  field,  equation  (3-62),  shows 
the  same  z-depencence  for  the  density  distribution. 

A  second  zone  in  this  plasma  is  that  in  the  vicinity  of  the  center 
of  the  field  coils  in  which  wall,  losses  due  to  radial  diffusion  to 
the  walls  control  the  density  distribution.  Here 

■H(r)  \  no  X(^r)  =  lio  X  (2~]=r  r)  ( 5-2) 

where  n0  is  the  density  in  the  center  of  the  cylinder. 

Then,  as  a  first  approximation,  we  shall  assume  that  the  £on£za~ 
tion  coefficient,  V    9   in  the  zone  of  parallel  diffusion  will  not  ciffer 
significantly  from  that,  of  the  corresponding  stable  radial  configura= 
tion.  In  higher  approximations  we  would  have  to  calculate  the  change 
in  this  quantity  due  to  induced  changes  it   electron  temperature  by 
the  different  fields,  but  we  shall  assume  these  to  be  negligible. 
Consequently, 

yVp  (5=3) 

Therefore,  we  obtain  from  equations  (5-1)  and  (5~2) 
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The  surface  separating  the  tvo  models  (equation  (3-lil))  is  that  at 
vh  i  ch 


„ 


i 


i' 


U  S^Cjd) 


D 


H\ 


(S-S) 


Substitution  from  equation  (5~U)  reduces  this  to 


'  Cct(f^) 


Jd^^Y 


(5-6) 


Figure  (5~?)  is  a  nomogram  from  which  the  z-coordinate  of  the 
separation  surface,  z$  can  be  obtained  as  a  fraction  of  z  ,  using 
r  (as  a  fraction  of  cylinder  radius  R)  and 


as  arguments.  From  equation  (5«=7)»  it  is  apparent  that  "a"  decreases 

(with  B)  along  the  cylinder  axis  (indicated  on  Figure  (5=1))  beyond 

the  end  of  the  field  coils.  Consequently,  when  using  the  nomogram 

to  determine  zQ,  it  is  necessary  to  use  the  value  of  "a"  at  the 

junction  of  the  tvo  models.  This  requires  a  trial-and<=error  procedure, 

In  order  to  illustrate  this  method  of  approximation,  and  to 

estimate  these  effects,  a  system  with  the  fol loving  dimensions  vas 

selected? 

Plasma  container  radius?    2.5  cm. 

magnetic  coils i  width  -  15  cm. 

spacing  -  5  cm. 

number  -  10 

inside  diameter  -  25  cm. 

outside  diameter  -  50  cm. 

overall  length  -  195  cm. 
10U 


The  axial  field  of  the  above  array  of  field  coils  vas  computed  with 
a  digital  computer  (figure  5-3) ,  using  the  thick  coil  equations  fran 
reference  [_20]    .     The  author  states  that     they  agree  to  vithin  a 
few  percent  of  measured  values,  and  for  the  region  vithin  the  coil 
array3   experimental  curves  from  a  very  similar  array  were  available 
and  agreed  very  well.     The  field  was  assumed  constant  through  any 
cross-section  of  the  cylinder,  which  is  a  good  approximation  because 
of  the  small  radius  of  the  cylinder.     Figure  5-1  shows  the  field  lines 
as  computed  from  the  field  strength  data. 

Figure  (5-1)  also  shovs  the  separation  of  the  two  zones  of 

2 
diffusion  for  three  different  orders  of  magnitude  of  B  .,   namely, 

ol    _     /OOO-I 

Tl^         10,000-1 

(5-8) 
It  is  to  be  noted  that  for  B  *  i£-   (a  «■  10)  the  effect  of  the 

/v- 

parallel  diffusion  is  very  small,  and  is  centered  at  the  end  of 

the  coil  array.  For  B  *=•  jj-f:    (a  »  31.6)  the  effect  is  appreciable^ 

and  extends  more  than  half-way  into  the  cylinder.  Figure  (5-k) 

is  a  density  profile  of  the  resulting  particle  distribution^  Figure 

(5 -5) >   selected  cross-sections  from  Figure  (5-1) I  and,  Figure  (5-6) 

is  a  plot  of  constant  density  curves. 

One  feature  of  the  density  distribution  vhich  must  be  discussed 

is  the  question  of  determining  the  density  distribution  of  the  central 

zone  when  it  does  not  extend  to  the  walls  but  rather  terminates  at 

the  edge  of  a  longitudinal  zone.  A  longitudinal  current  (between 
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electrodes  at  either  end  of  the  cylinder)  is  responsible  for  the 
ionization  within  the  plasma,  and  this  current  must  be  constant 
through  the  plasma.  Further,  we  have  assumed  that  the  ionization 
coefficient,  p  ,  is  constant.  From  this  it  follows  that  the  electron 
temperature  is  constant,  and  therefore  the  longitudinal  electric  field 
must  be  approximately  constant.  But  if  the  field  and  current  are 
approximately  constant,  then  the  total  number  of  particles  in  a 
cross-section  must  be  constant.  Further,  the  plasma  in  the  radial- 
diffusion  region  vhich  is  not  controlled  by  the  cylinder  boundaries 
because  of  the  intervening  region  of  parallel  diffusion  is  expected 
to  maintain  its  Bessel  distribution  (thus  satisfying  the  diffusion 
equation  and  not  requiring  a  different  ionization  coefficient).  As 
a  consequence  of  the  constant  field  and  current  restrictions,  however, 
the  total  cross-section  particle  density  (particles  per  unit  length^ 

is  a  constant,  equal  to  the  total  number  of  particles  per  unit  length 
at  the  center  of  the  field  coil.  Therefore  the  normalization  of  the 
Bessel  distribution  for  the  center  of  the  plasma  in  regions  where  the 
distribution  at  greater  radii  is  controlled  by  longitudinal  diffusion 
must  be  so  determined  that  equation  (5-9)  is  satisfied. 

One  effect  vhich  we  ignore  in  the  previous  formulation  is  that 
due  to  the  small  radial  component  of  the  magnetic  field  at  the  end  of  the 
coil.  This  component  is  strongest  at  greater  radii,  and  decreases 
the  mobility  of  the  particles  in  the  longitudinal  direction.  Con- 
sequently, it  would  contribute  further  to  the  constriction  of  the 
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plasma  at  this  point.  VJe  neglect  this  influence,  assuming  the 
curvature  of  the  field  lines  to  be  very  anall  (recall  that  the  radial 
scale  of  Figure  (5-1)  is  greatly  enlarged). 

It  should  be  noted  that  there  is  a  discontinuity  in  the  direction 
of  V  n  along  all  the  dividing  surfaces  between  zones,  and  a  discon° 
tinuity  in  particle  density,  n,  along  the  a  ■  31.6  boundary  from  90  cm 
to  130  cm  and  the  a  -  100  boundary  from  the  center  to  l6h  cm.  The 
former  discontinuity  is  a  consequence  of  our  approximation.  The 
density  discontinuity  occurs  because  the  dividing  surface  lies  along 
one  of  the  field  lines,  and  is  a  result  of  the  fact  that  the  sine  dis=» 
tribution  for  the  longitudinal  model  cannot  be  made  to  match  the 
normalized  radial  distribution  along  such  a  line  except  at  the  end 
which  is  nearer  the  end  of  the  coil.  The  transition  regions  are 
centered  along  these  dividing  surfaces,  and  are  the  region  of  greatest 
uncertainty.  Along  those  surfaces  at  which  the  particle  density,  n, 
is  discontinuous,  this  approximation  is  even  more  uncertain.  Because 
of  this,  the  a  ■  100  description,  with  a  discontinuity  in  both 
and  n  along  its  entire  length,  must  be  much  too  crude  to  be  valid, 
and  further  attempts  to  fit  this  approximation  to  the  effects  of 
such  large  fields  are  not  attempted.  The  description  for  a  ■  31.6 
is  believed  to  be  a  valid  "order-of-magnitude1*  approximation. 

In  Figure  (5-7)  the  heat  liberated  at  the  walls  by  the  ions  and 

electrons  carried  to  them  is  indicated  as  a  function  of  position  along 

the  cylinder.  There  are  three  curves;  all  calculated  for  a  ■»  31.6s 

(1)  the  wall  heat  loss  for  normal  ambipolar  diffusion  without  a 

magnetic  field,  (2)  the  •'all  heat  loss  for  radial  diffusion  with  a 

magnetic  field,  accounting  for  the  increase  in  radial  diffusion  vhich 

accompanies  the  decrease  in  magnetic  field  outside  the  coils,  and 
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(3)  the  estimated  heat  loss  for  the  magnetic  field,  taking  into  account 
the  longitudinal  as  well  as  radial  diffusion.  This  latter  curve  is  a 
very  crude  approximation  based  on  the  results  predicted  by  this  model 
with  adjustments  to  compensate  for  the  errors  introduced  by  the 
discontinuities  at  the  separation  surfaces.  The  figure  shows  that 
the  heating  at  the  wall  vrithin  the  region  of  parallel  diffusion  and 
within  the  coil  array  is  greatly  decreased  because  the  particles  are 
transported  to  the  end  of  the  cylinder  instead  of  directly  to  the  walls. 
However,  outside  the  coil,  where  particles  from  within  the  coil  are 
carried  by  the  longitudinal  diffusion,  the  effects  of  such  additional 
heating  is  masked  by  the  heating  due  to  the  increased  raidal  current 
which  flows  because  of  the  lower  magnetic  field  intensity. 

The  values  of  magnetic  field  B  for  vfaich  this  effect  can  be 
approximately  described  by  this  model  are  those  whose  order  of  magni<= 
tude  corresponds  to  a  «*31.6.  Smaller  fields  (a  ^10)  will  not  cause 
an  appreciable  end-effect  constriction.  Larger  fields  (  (X   2^  100) 
are  not  describabJe  by  this  model.  For  a  plasma  which  fits  the  sample 
data  given  in  section  3-B-l,  this  is  seen  from  equation  (5-8)  to  be 
approximately 


Q  *  JTEP      z*     L\>         weberA2   (5-10) 


It, 
It  must  be  recalled,  however,  that  the  above  description  is  based 

upon  the  walls  being  perfect  insulators.  Electron  diffusion,  pre- 
dominately parallel  to  the  field  lines  because  of  the  small  electron 
Larmor  radius,  is  restrained  by  the  build-up  of  a  negative  wall  charge 
in  the  vicinity  of  the  magnetic  field  line  intersections  with  the 
cylinder  walls.  Similarly,  the  ion  diffusion,  which  would  be  predominantly 
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radial  because  of  the  larger  Larmor  radii,  would  be  opposed  by  a  large 
positive  wall  charge  within  the  coil.  In  this  manner,  the  conditions 
of  equal  electron  and  ion  radial  current  and  equal  longitudinal  current 
are  met.  However,  it  is  known  that  insulators  have  a  surface  conduc- 
tivity which  depends  on  the  material  and  the  temperature.  The  effect 
of  this  conductivity  is  to  permit  an  exchange  of  carriers  along 
the  wall.  Consequently,  we  obtain  diffusion  which  is  somewhere  between 
Allis  diffusion  [l%]   and  our  qua si-ambi polar  diffusion.  The  wall 
charge  build-up  is  avoided  by  the  exchange  of  ions  and  electrons 
between  negatively  and  positively  charged  zones.  This  exchange  results 
in  recombination  along  the  length  of  the  cylinder  and  would  probably 
invalidate  the  conclusions  of  Figure  (5-7)  if  the  currents  were 
appreciable.  Further,  if  this  surface  conductivity  were  sufficiently 
high,  the  diffusion  would  approach  Allis  diffusion,  under  which  condi- 
tion the  ratio  of  the  parallel  and  perpendicular  components  of  the 
diffusion  tensor  would  no  longer  be  1  +  /4^/4-b  v     but  would  be 
1  +  /. L  B,     Thus  the  field  strength  which  would  produce  a  con- 
striction describable  by  this  model  would  be  reduced  to 

The  description  of  the  density  distributions  would  not  be  affected 
by  the  change  in  diffusion  process  since  the  only  effect  on  the 
diffusion  is  a  change  in  the  magnitudes  of  the  components  of  the 
diffusion  tensor.  In  practice  the  diffusion  would  be  somewhere  be- 
tween the  extremes  of  no  surface  conductivity  and  perfect  conductivity. 
Therefore,  the  value  of  B  to  which  these  calculations  correspond  would 
intermediate  between  OJt  weber/m^  and  h  weber/m^,  depending  upon 

the  surface  current  conducted  by  the  plasma  container  walls. 
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It  should  be  noted  that  this  description  predicts  negligible  end  effects 
for  the  plasma  described  for  fields  below  about  0.1  weber/m  .  For 
particles  vith  higher  mobilities,  however,  this  effect  could  occur  in 
smaller  magnetic  fields. 
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APPENDIX   I  -  Curvilinear  Coordinates?  Transformation  Equations 

Because  of  notation  differences  between  references,  the  trans- 
formation equations  for  all  coordinate  systems  employed  in  this  thesis 
are  presented  belows 

1.  Cylindrical  Coordinates  (r,tf,z)   -  Figure  1-1 

x  =  r  cos^ 

y  =  r  sinflf  (l=l) 

z  =  z 

'•-here  r  |.  0,     O£0fi2TT  -oo<z<   oo 

2.  Spherical  Coordinates     (r,9,(2Q    -  Figure  1-2 

x  =  r  sin9  cos$ 

y  =  r  sir9  cos^f                            (l™?) 

z  ~   r  cos9 

where  r  >  0  0  £  $  :£  2  "H    0  £  9  ^  "FT 

3.  Prolate  Spheroidal  Coordinates  (71 ,  ^,  y)     «=  Figure  1=3 


*f  ? 


where  -1  1  tf  <  1    1    ?^ao     0  <.#  £  2tj 
In  the  prolate  spheroidal  system  the  surface  *T  «=  constant  >1  is  an 
elongated  ellipsoid  of  revolution  with  major  axis  of  length  d?"  and 
minor  axis  of  length  d  (  ^  -l)2.  The  degenerate  surface  Z   s  1  is  the 
straight  line  along  the  z-axis  from  z  =  -jd  to  z  =  +^do  The  surface 
|V/|  =  constant  <1  is  a  hyperboloid  of  revolution  of  two  sheets  with  an 
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asymptotic  cone  vhose  generating  line  passes  through  the  origin  and 
is  inclined  at  the  '/3  ~   cos"1/!   to  the  z-axis.  The  degenerate  surface 
I'M  I    =  1  is  that  part  of  the  z-axis  for  which  jz|  <  jd.  The  surface 
g!  =  constant  is  a  plane  through  the  z-axis  forming  the  angle  4  vith  the 
x,z-plane. 

It.  Oblate  Spheroidal  Coordinates  (jf,  J  ,$)     -  Figure  I=\\ 

I   vhere  -1  *  VI    <.  1    0  £.  ^  <£.  co     0  i=  #  £•  21T 

In  the  oblate  spheroidal  system  the  surface  m  ■  constant  >  0  is  a 

+1)2 
and  minor  axis  of  length  d|F|  •  The  surface  "?  =0  is  a  circular  disk 
of  radius  a  =  -|-d  which  lies  in  the  x,y-plane  and  is  centered  at  the 
origin.  The  surface  \y\\-   constant  ^  1  is  a  hyperboloid  of  revolution 
of  one  sheet  with  an  asymptotic  cone  whose  generating  line  passes  through 
the  origin  and  is  inclined  at  the  angle  &   ■  cos™  /I    to  the  z=axisa 
The  degenerate  surface  N~l  is  the  z-axis.  The  surface  Yl  K)  is  the 
x,y-plane  except  for  the  circular  disk  T=0.  The  surface  $  -   constant 
is  again  the  plane  through  the  z-axis  taking  the  angle  $  vith  the  x, 
z -plane.  The  transformation  from  prolate  to  oblate  coordinates  is 
made  by  replacing  V    by  iK     in  the  prolate  equations0 

Reference  for  spheroidal  systems?  Flammer,  Spheroidal  Wave  Functions  6 
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figure  X-l      Cuh^nAncA}   Ccord\r\oM&         ftmre   J-Z    SpVierico-l    Goorijml&s 


figure  1-3     FV*»Ute    Spheroid*) 

Coordmtdes    -  axis  tff 


££,<Dr<Am<*l*e£    —axis    <sf 
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APPENDIX  II 
SPHEROIDAL  VAV£  FUNCTIONS 

The  particle  density  distribution  in  the  vicinity  of  a  prolate 
spheroid  is  describable  by  the  prolate  angle  function  of  the  first 

Sf  I)    f  V 

k,  (£.71)  and  the  prolate  radial  functions  of  the  first  and 

second  kinds,  Roi    ('/ f  )      and         Roi  L^A)        >  respectively0 

Flammer  [6\   has  tabulated    Soj  (C,?f)  (See  Fi9»  II-l,  which  is 

taken  from  his  data)  and  the  radial  functions  Roi(C,£)   and  j^oi  C^,^) 

for  different  values  of  c  and  the  independent  variables  >'l   and  £  . 

The  range  of  T     for  the  radial  functions  in  his  tables  is,  however, 

quite  restricted  (  <    =1.0C£,  1.020,  l.Oltli,  and  1.077)  and  consequently 

it  vas  necessary  to  compute  further  tables  of  the  radial  functions 

(Table  II-l  and  II-2). 

To  compute      No/C^y     and      f\oi(y\)the  following  expansions 


in  terms  of  spherical  Bessel  functions  and  Legendre  Polynomials  vere 
obtained  from  Flammer's  general  equations.  (lnl.l£)  for  ^rr\v\(c/\)and 
(U.2.6)for     f?$fc,^ 

,0)         v  _r '  .r-t  fit. 

ox   (C,%)    =    Y 

and  r*  ' 


(?!'.  fc  0  -  Yr~-Lur<tr^  1rK)      (n-D 


£D 


where  from  Flammer1  s  equation  (U.  1.8)  and  (U.2.£b)s 


and 


and  where  the  prime  over  the  summation  sign  indicates  that  the 

summation  is  over  only  odd  values  of  r,  and  vhere   n  <-  (c)  and  c\pj,(c) 

121* 


are  e>pansion  coefficients.  Flammer  tabulates  these  coefficients 
to  five  significant  figures. 

The  prolate  radial  functions  of  the  first  kinds   ^0(  (c,^)  (Table 
II-l),  were  calculated  using  a  desk  calculator  and  the  Mathematical 
Tables  Project  Tables  of  Spherical  Bessel  Functions  JjLi]  „  The  figures 
have  been  spot  checked,  and  are  believed  correct  to  plus  or  minus  one 
in  the  last  digit. 

The  Prolate  Radial  Functions  of  the  second  kind,  ^(C5)(Table  II-2) 
were  calculated  on  the  CI£»l60li  Computer  using  equation  £ll=>2\  Most 
of  the  Legendre  Functions  required  for  the  expansions  were  also  genera- 
ted on  the  computer,  using  the  following  equations  given  by  Jahnke 
and  Emde  QLO]  . 

a   tc\     l'3S'"(2w'l)^  J^^.JMNM  ^   -.   x 

ft*  (?)  s  —^7 (f  -5^  *J*&jfc-$  -••■1(ji-3) 

Because  of  the  slow  convergence  of  the  equations  (II -3)  for  large 
values  of  £"  only  those  functions  for   X  <   ),2  were  computed 0  The 
values  for   X    >     1.2  were  taken  from  Tables  of  Associated  Legendre 
Functions  (J.  2]  .  The  Values  of  Q0  for  7  <  1.2  were  calculated  on 
a  desk  calculator  using  the  following  relation  from  Jahnke  and  Emde  [103« 

Qe(f)  =  esKTYf)  =T&4^f  n-5 

and  tables  of  logarithms. 

The  radial  functions   j^ol(c^)and  ]*\6|  (c;^)  are  tabulated  here 

for  values  of  c  between  0.1  and  2.0  and  for  values  of   T  between  1.0 

) 

and3.6  and  1.0  and  2.0  respectively.  These  functions  are  graphed  in 

Figures  (II-2),  (II-3)«>  and  (II-lt).  The  tabulated  values  are  correct 
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to  the  number  of  digits  shown,  except  that  all  five  digit 
results  are  considered  correct  to  plus  or  minus  one  unit  in 
the  last  digit  since  Flammer's  coefficients  are  listed  to  only- 
five  significant  digits.  The  decreasing  number  of  significant 
digits  for  higher  values  of  c  and  ^  is  a  consequence  of  the 
decreasing  rate  of  convergence  of  the  series  (II-2)  for  rela- 
tively large  c  and  J  .  The  number  of  significant  digits  was 
determined  by  inspection  of  the  terms  and  generous  estimates  cf 
the  maximum  magnitude  of  the  neglected  terms.  The  author  has 
endeavered  to  err  on  the  conservative  side  in  reporting  these 
values,  and  believes  all  digits  reported  to  be  correct,. 
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